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The very special characteristic of the proximity effect in superconductor-ferromagnet systems is 
the damped oscillatory behavior of the Cooper pair wave function in a ferromagnet. In some sense, 
this is analogous to the inhomogeneous superconductivity, predicted long time ago by Larkin and 
Ovchinnikov (1964), and Fulde and Ferrell (1964), and constantly searched since that. After the 
qualitative analysis of the peculiarities of the proximity effect in the presence of the exchange field, 
the author provides a unified description of the properties of the superconductor-ferromagnet het- 
erostructures. Special attention is paid to the striking non-monotonous dependance of the critical 
temperature of the multilayers and bilayers on the ferromagnetic layer thickness and conditions 
of the realization of the 'V- Josephson junctions. The recent progress in the preparation of 
the high quality hybrid systems permitted to observe on experiments many interesting effects, 
which are also discussed in the article. Finally, the author analyzes the phenomenon of the 
domain-wall superconductivity and the influence of superconductivity on the magnetic structure 
in superconductor-ferromagnet bilayers. 
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I. INTRODUCTION 

Due to their antagonistic characters, singlet superconductivity and ferromagnetic order cannot coexist in bulk 
samples with realistic physical parameters. Ginzburg (1956) was the first to set up theoretically the problem of 
magnetism and superconductivity coexistence taking into account the orbital mechanism of superconductivity de- 
struction (interaction of the superconducting order parameter with a vector-potential A of the magnetic field). After 
the creation of BCS theory, it became clear that superconductivity (in the singlet state) can be also destroyed by the 
exchange mechanism. The exchange field, in the magnetically ordered state, tends to align spins of Cooper pairs in the 
same direction, thus preventing a pairing effect. This is the so-called paramagnetic effect (Saint- James et al., 1969). 
Anderson and Suhl (1959) demonstrated that ferromagnetic ordering is unlikely to appear in the superconducting 
phase. The main reason for that is the suppression of the zero wave- vector component of the electronic paramagnetic 
susceptibility in the presence of superconductivity. In such situation the gain of energy for the ferromagnetic ordering 
decreases and instead of the ferromagnetic order the non-uniform magnetic ordering should appear. Anderson and 
Suhl (1959) called this state cryptoferromagnetic. 

The 1977 discovery of ternary rare earth (RE) compounds (RE)Rli4B4 and (RE)Mo6Xg (X=S, Se) (as a review see, 
for example, Maple and Fisher, 1982) provided the first experimental evidence of magnetism and superconductivity 
coexistence in stoichiometrical compounds. It turned out that in many of these systems, superconductivity (with 
the critical temperature T c ) coexists rather easily with antiferromagnetic order (with the Neel temperature T/v), and 
usually the situation with Tjv < T c is realized. 

The more recent discovery of superconductivity in the quaternary intermetallic compounds (RE)Ni 2 B 2 C (as a review 
see, for example, Miiller and Narozhnyi, 2001) gives another example of antiferromagnetism and superconductivity 
coexistence. 

Indeed, superconductivity and antiferromagnetism can coexist quite peacefully because, on average, at distances of 
the order of the Cooper pair size (superconducting coherence length) the exchange and orbital fields arc zero. Much 
more interesting a re-entrant behavior of the superconductivity was observed in ErRh4B4 and HoMogSs (Maple 
and Fisher, 1982). For example, ErRh 4 B4 becomes superconductor below T c = 8.7 K. When it is cooled to the 
Curie temperature « 0.8 K an inhomogeneous magnetic order appears in the superconducting state. With further 
cooling the superconductivity is destroyed by the onset of a first-order ferromagnetic transition at the second critical 
temperature T c2 ~ 0.7 K. HoMo 6 Ss gives another example of the re-entrant superconductivity with T c = 1.8 K, 
6 w 0.74 K, and T c2 w 0.7 K. 

In these compounds at Curie temperature, following the prediction of Anderson and Suhl (1959) a non-uniform 
magnetic order appears. Its presence was confirmed by neutron scattering experiments. The period of this magnetic 
structure is smaller than the superconducting coherence length, but larger than the interatomic distance. In some 
sense this structure is a realization of the compromise between superconductivity and ferromagnctism : for the super- 
conductivity it is seen as an antiferromagnetism, but for the magnetism it looks like a ferromagnetism. Theoretical 
analysis, taking into account both orbital and exchange mechanisms and magnetic anisotropy (as a review see Bu- 
laevskii et al., 1985), revealed that the coexistence phase is a domain- like structure with very small period. The region 
of magnetism and superconductivity coexistence in ErRli4B4 and HoMoeSg is narrow, but in HoMogSeg the domain 
coexistence phase survives till T = K. 

The first truly ferromagnetic superconductors £/Ge 2 (Saxena et al., 2000) and URhGe (Aoki et al., 2001) have been 
discovered only recently, and apparently the coexistence of superconductivity with ferromagnetism is possible due to 
the triplet character of the superconducting pairing. Indeed, the superconductivity in URhGe (Aoki et al., 2001) 
appears below 0.3 K in the ferromagnetic phase which has the Curie temperature = 9.5 K; this makes the singlet 
scenario of superconductivity rather improbable. 

Though the coexistence of singlet superconductivity with ferromagnctism is very unlikely in bulk compounds, it may 
be easily achieved in artificially fabricated layered ferromagnet /superconductors (F/S) systems. Due to the proximity 
effect, the Cooper pairs can penetrate into the F layer and induce superconductivity there. In such case we have 
the unique possibility to study the properties of superconducting electrons under the influence of a huge exchange 
field acting on the electron spins. In addition, it is possible to study the interplay between superconductivity and 
magnetism in a controlled manner, since varying the layer thicknesses we change the relative strength of two competing 
orderings. The behavior of the superconducting condensate under these conditions is quite peculiar. 

Long time ago Larkin and Ovchinnikov (1964), and Fulde and Ferrell (1964) demonstrated that in a pure ferro- 
magnetic superconductor at low temperature the superconductivity may be non-uniform. Due to the incompatibility 
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of ferromagnetism and superconductivity it is not easy to verify this prediction on experiment. It occurs that in S/F 
systems there exists some analogy with the non-uniform superconducting state. The Cooper pair wave function has 
damped oscillatory behavior in a ferromagnet in contact with a superconductor. It results in many new effects that 
we discuss in this article : the spacial oscillations of the electron's density of states, the non-monotonous dependance 
of the critical temperature of S/F multilayers and bilayers on the ferromagnet layer thickness, the realization of 
the Josephson "71""- junctions in S/F/S systems. The spin-waive effect in the complex S/F structures gives another 
example of the interesting interplay between magnetism and superconductivity, promising for the potential appli- 
cations. We discuss also the issues of the localized domain- wall superconductivity in S/F bilayers and the inverse 
influence of superconductivity on ferromagnetism, which favors the non-uniform magnetic structures. An interesting 
example of atomic thickness S/F multilayers is provided by the layered superconductors like Smi.ssCeo.isCuC^ and 
RuSr 2 GdCu 2 08. For such systems the exchange field in F layer also favors the 'V-phase behavior, with an alternating 
order parameter in adjacent superconducting layers. 

Note that practically all interesting effects related with the interplay between the superconductivity and the mag- 
netism in S/F structures occurs at the nanoscopic range of layers thicknesses. The observation of these effects became 
possible only recently due to the great progress in the preparation of high-quality hybrid F/S systems. The exper- 
imental progress and the possibility of potential applications in its turn stimulated a revival of the interest to the 
superconductivity and ferromagnetism interplay in heterostructures. It seems to be timely to review the present state 
of the research in this domain and outline the perspectives. 



II. PARAMAGNETIC LIMIT AND QUALITATIVE EXPLANATION OF THE NON-UNIFORM PHASE FORMATION 
A. The (H, T) phase diagram 

For a pure paramagnetic effect, the critical field of a superconductor H p at T = may be found from the comparison 
of the energy gain AE n due to the electron spin polarization in the normal state and the superconducting condensation 
energy AE S . Really, in the normal state, the polarization of the electron gas changes its energy in the magnetic field 

by 

A£„ = -Xny, (1) 

where x« = 2^^iV(0) is the spin susceptibility of the normal metal, \ib is the Bohr magneton, 27V(0) is the density 
of electron states at Fermi level (per two spin projections), and the electron g factor is supposed to be equal to 2. 
On the other hand, in a superconductor the polarization is absent, but the BCS pairing decreases its energy by 

AE S = -N(0)^, (2) 

where A = 1.76T C is the superconducting gap at T = 0. From the condition AE n — AE S , we find the Chandrasckhar 
(1962) - Clogston (1962) limit (the paramagnetic limit at T = 0) 

H p (0) = -^-. (3) 

Note that it is the field of the first-order phase transition from a normal to a superconducting state. The complete 
analysis (Saint-James et at, 1969) demonstrates that at T = this critical field is higher than the field of the second 
order phase transition H^ 1 (0) = Ao/2/zb, and the transition from a normal to a uniform superconducting state is of 
the second-order at T* < T < T c only, where T* = 0.56T C , H* = H{T*) = 0.61A / >b = 1.05T c //x B . However, Larkin 
and Ovchinnikov (1964), and Fulde and Ferrell (1964) predicted in the framework of the model of pure paramagnetic 
effect the appearance of the non-uniform superconducting state with a sinusoidal modulation of the superconducting 
order parameter at the scale of the superconducting coherence length £ s (the FFLO state). In this FFLO state, the 
Cooper pairs have a finite momentum, compared with zero momentum in conventional superconductors. Recently 
Casalbuoni and Nardulli (2004) reviewed the theory of the inhomogeneous superconductivity applied to the condensed 
matter and quantum chromodynamics at high density and low temperature. 

The critical field of the second-order transition into FFLO state goes somewhere above the first-order transition line 
into a uniform superconducting state (Saint-James et at, 1969). At T = 0, it is H FFLO (0) = 0.755A /^b (whereas 
Hp = 0.7Aq//xb). This FFLO state only appears in the temperature interval < T < T* , and is sensitive to impurities 
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(Aslamazov, 1968). In a dirty limit it is suppressed, and the first-order transition into the uniform superconducting 
state takes place instead. The phase diagram for the 3D superconductors in the model of pure paramagnetic effect 
is presented in Fig. 1 (Saint- James et al, 1969). Up to now, there were no unambiguous experimental proofs of this 
state observation. Note however that, recently, the magnetic-field-induccd superconductivity has been observed in 
the quasi two-dimensional organic conductor {BETS)-2FeCli (Uji et al, 2001) which is an excellent candidate for 
the FFLO state formation (Balicas et al, 2001 and Houzet et al, 2002). 



B. Exchange field in the ferromagnet 



In a ferromagnet the exchange interaction between the electrons and the magnetic moments may be considered as 
some effective Zeeman field. In the case of magnetic moments with spin Si, localized in the sites r i; their interaction 
with electron spins is described by the exchange Hamiltonian 



= J d 3 r*+(r) J(r - r^cr j * 



H mt = / d 3 r*+(r) <yj(r- r^a *(r), (4) 



where ^(r) is the electron's spinor operator, a = {a x ,a y , a z } are the Pauli matrices, and J(r) is the exchange integral. 
Below the Curie temperature 6, the average value of the localized spins (Si) is non-zero, and the exchange interaction 
may be considered as some effective Zeeman field H e ^ = J J(r)d 3 r, where n is the concentration of localized 

moments, and the spin quantization z-axis is chosen along trie ferromagnetic moment. It is convenient to introduce 
the exchange field h as 

h = f-i B H eff = (5?) n j J(r)d 3 r = s(T)h , (5) 

where s(T) = (Sf) / (S?) T _ is the dimensionless magnetization and ho is the maximum value of an exchange field at 
T = . The exchange field h describes the spin-dependent part of the electron's energy and the exchange Hamiltonian 
Eq. (4) is then simply written as 

H int = I d 3 r^+{v)ha^(v). (6) 



If we also want to take into account the proper Zeeman field of magnetization M, then we may simply replace h in 
Eq. (6) by h +4irMng . The reader is warned that in principle, if the exchange integral is negative, the exchange 
field may have the direction opposite to the magnetic moments and the interesting compensation Jaccarino-Peter 
(1962) effect is possible. However, in the ferromagnetic metals, the contribution of the magnetic induction to the spin 
splitting is several order of magnitude smaller than that of the exchange interaction and may be neglected. In the case 
of the Rudcrman-Kittcl-Kasuya-Yosida (RKKY) mechanism of the ferromagnetic ordering, the Curie temperature 
6 <~ Hq/Ef and in all real systems the exchange field ho » 6, T c . This explains that the conditions of singlet 
superconductivity and ferromagnetism coexistence are very stringent. Indeed, if 6 > T c the exchange field in a 
ferromagnet h » T c , which strongly exceeds the paramagnetic limit. On the other hand, if 6 < T c then, instead of 
the ferromagnetic transition the inhomogeneous magnetic ordering appears (Maple and Fisher, 1982; Bulaevskii et al, 
1985). The very high value of the exchange field in ferromagnet permits us to concentrate on the paramagnetic effect 
and neglect the orbital one (note that well below the Curie temperature the magnetic induction A-kM in ferromagnets 
is of the order of several koe only). 



C. Why does the Fulde-Ferrell-Larkin-Ovchinnikov state appear? 

What is the physical origin of the superconducting order parameter modulation in the FFLO state ? The appearance 
of modulation of the superconducting order parameter is related to the Zeeman's splitting of the electron's level under 
a magnetic field acting on electron spins. To demonstrate this, we consider the simplest case of the ID superconductor. 

In the absence of the field, a Cooper pair is formed by two electrons with opposite momenta +kp and — kp 
and opposite spins ( j) and (J.) respectively. The resulting momentum of the Cooper pair fcp + {—kp) = 0. Under a 
magnetic field, because of the Zeeman's splitting, the Fermi momentum of the electron with spin (f ) will shift from kp 
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to k\ = kp + Skp , where Skp = jisH/vp and vp is the Fermi velocity. Similarly, the Fermi momentum of an electron 
with spin (J.) will shift from — kp to k 2 = —kp + Skp (see Fig. 2) . Then, the resulting momentum of the Cooper 
pair will be k\ + fc 2 = 2Skp ^ 0, which just implies the space modulation of the superconducting order parameter 
with a resulting wave- vector 25kp. Such type of reasoning explains the origin of the non- uniform superconducting 
state formation in the presence of the field acting on electron spins, and, at the same time, demonstrates the absence 
of a paramagnetic limit (at 0) for the ID superconductor (Buzdin and Polonskii, 1987). For 3D (Larkin and 
Ovchinnikov, 1964 and Fulde and Ferrcll, 1964) or 2D (Bulaevskii, 1973) superconductors, it is not possible to choose 
the single wave vector Skp which compensates the Zeeman splitting for all electrons on the Fermi surface (as Skp 
depends on direction of vp), and the paramagnetic limit is preserved. However, the critical field for a non- uniform 
state at T = is always higher than for a uniform one. However, the critical field for a non-uniform state at T = 
is always higher than for a uniform one. At finite temperature (when T > /j-bH ), the smearing of the electrons 
distribution function near the Fermi energy decreases the difference of energies between the non-uniform and uniform 
states. As it follows from the microscopical calculations, at T > T* = 0.56T C the uniform superconducting phase is 
always more favorable (Saint-James et al, 1969). 



D. Generalized Ginzburg Landau functional 

Qualitatively, the phenomenon of the FFLO phase formation and the particularities of the proximity effect in S/F 
systems may be described in the framework of the generalized Ginzburg-Landau expansion. Let us first recall the 
form of the standard Ginzburg-Landau functional (see, for example, De Gennes, 1966) 



F = a|^| 2 +7 



(7) 



where tp is the superconducting order parameter, and the coefficient a vanishes at the transition temperature T c . At 
T < T c , the coefficient a is negative and the minimum of F in Eq. (7) is achieved for a uniform superconducting state 
with |V! = -f • If 

we consider also the paramagnetic effect of the magnetic field, all the coefficients in Eq. (7) will 
depend on the energy of the Zeeman splitting /j-bH , i. e. an exchange field h in the ferromagnet. Note that we neglect 
the orbital effect, so there is no vector-potential A in Eq. (7). To take into account the orbital effect in the Ginzburg- 
Landau functional, we may substitute the gradient by its gauge-invariant form y — > y — ^A. Usually, the orbital 
effect is much more important for the superconductivity destruction than the paramagnetic one. It explains why in the 
standard Ginzburg-Landau theory there is no need to take into account the field and temperature dependence of the 
coefficients 7 and b. However, when the paramagnetic effect becomes predominant, this approximation fails. What 
are the consequences ? If it was simply some renormalization of the coefficients in Ginzburg-Landau functional, the 
general superconducting properties of the system would basically be the same. However, the qualitatively new physics 
emerges due to the fact that the coefficient 7 changes its sign at the point (H* ,T*) of the phase diagram, see Fig. 1. 
The negative sign of 7 means that the minimum of the functional does not correspond to an uniform state anymore, 
and a spatial variation of the order parameter decreases the energy of the system. To describe such a situation it is 
necessary to add a higher order derivative term in the expansion (7), and the generalized Ginzburg-Landau expansion 
will be: 



F G =a(H,T)\1>\ 2 +7(H,T) 



+ 



ri(H,T) 



b(H,T) 



+ 



(8) 



The critical temperature of the second order phase transition into a superconducting state may be found from the 
solution of the linear equation for the superconducting order parameter 



aip-jAip + ^A 2 V> = 0. 



(9) 



If we seek for a non-uniform solution ip 
vector q and is given by the expression 



ip exp(iqr), the corresponding critical temperature depends on the wave- 
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V 4 
2 q - 



(10) 
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Note that the coefficient a may be written as a — a(T — T CU (H)), where T CU (H) is the critical temperature of the 
transition into the uniform superconducting state. In a standard situation, the gradient term in the Ginzburg-Landau 
functional is positive, 7 > 0, and the highest transition temperature coincides with T CU (H); it is realized for the 
uniform state with q = 0. However, in the case 7 < 0, the maximum critical temperature corresponds to the finite 
value of the modulation vector q^ = —7/^ and the corresponding transition temperature into the non-uniform FFLO 
state T ci (H) is given by 

a = a(T cl -T cu ) = 1-. (11) 
Zrj 

It is higher than the critical temperature T cu of the uniform state. Therefore, we see that the FFLO state appearance 
may simply be interpreted as a change of the sign of the gradient term in the Ginzburg-Landau functional. A more 
detailed analysis of the FFLO state in the framework of the generalized Ginzburg-Landau functional shows that it 
is not an exponential but a one dimensional sinusoidal modulation of the order parameter which gives the minimum 
energy (Buzdin and Kachkachi, 1997; Houzet et at, 1999). In fact, the generalized Ginzburg-Landau functional 
describes new type of superconductors with very different properties, and the whole theory of superconductivity must 
be redone on the basis of this functional. The orbital effect in the framework of the generalized Ginzburg-Landau 
functional may be introduced by the usual gauge-invariant procedure V — > V — —A. The resulting expression for 
the superconducting current is quite a special one and the critical field may correspond to the higher Landau level 
solutions as well as new types of vortex lattices may exist (Houzet and Buzdin, 2000; Houzet and Buzdin, 2001). 



III. PROXIMITY EFFECT IN FERROMAGNETS 

A. Some generalities about superconducting proximity efFect 

The contact of materials with different long-range ordering modifies their properties near the interface. In the case 
of a superconductor-normal metal interface, the Cooper pairs can penetrate the normal metal at some distance. If 
the electrons motion is diffusive, this distance is of the order of the thermal diffusion length scale Lt ~ \J D /T. where 
D is the diffusion constant. In the case of pure normal metal the coresponding characteristic distance is £t ~ vf/T. 
Therefore the superconducting-like properties may be induced in the normal metal, and usually this phenomenon is 
called the proximity effect. At the same time the leakage of the Coopers pairs weakens the superconductivity near 
the interface with a normal metal. Sometime this effect is called the "inverse proximity effect", and it results in the 
decrease of the superconducting transition temperature in thin superconducting layer in contact with a normal metal. 
If the thickness of a superconducting layer is smaller than some critical one, the proximity effect totally suppresses the 
superconducting transition. All these phenomena and the earlier experimental and theoretical works on the proximity 
effect were reviewed by Deutscher and de Gennes (1969). 

Note that the proximity effect is a rather general phenomenon not limited by the superconducting phase transition. 
For example, in the case of the surface magnetism (White and Geballe, 1979) the critical temperature at the surface 
can be higher then the bulk one. In the result the magnetic transition at the surface induces the magnetisation nearby. 
On the other hand, the volume strongly affects the surface transition characteristics. 

However, the unique and very important characteristic of the superconducting proximity effect is the Andreev 
reflection revealed at the microscopical level. Andreev (1964) demonstrated how the single electron states of the 
normal metal are converted into Cooper pairs and explained the mechanism of the transformation at the interface 
of the dissipativc electrical current into the dissipationlcss supercurrent. An electron with an energy below the 
superconducting gap is reflected at the interface as a hole. The corresponding charge 2e is transferred to the Cooper 
pair which appears on the superconducting side of the interface. The manifestation of this double charge transfer is 
that for a perfect contact the sub-gap conductance occurs to be twice the normal state conductance. The classical 
work by Blonder, Tinkham and Klapwijk (1982) gives the detailed theory of this phenomenon. 

Andreev reflection plays a primary role for the understanding of quantum transport properties of superconduc- 
tor/normal metal systems. The interplay between Andreev reflection and proximity effect was reviewed by Pannetier 
and Courtois (2000). The reader can find a detailed description of the Andreev reflection in the normal metal- 
superconductor junctions in the framework of the scattering theory formalism in the review by Beenakker (1997). 
Recent review by Deutscher (2005) is devoted to the Andreev reflection spectroscopy of the superconductors. 
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B. Damped oscillatory dependence of the Cooper pair wave function in the ferromagnets 



The physics of the oscillating Cooper pair wave function in a ferromagnet is similar to the physics of the supercon- 
ducting order parameter modulation in the FFLO state - see section II. C. Qualitative picture of this effect has been 
well presented by Demler, Arnold, and Bcaslcy (1997). When a superconductor is in a contact with a normal metal 
the Cooper pairs penetrate across the interface at some distance inside the metal. A Cooper pair in a superconductor 
comprises two electrons with opposite spins and momenta. In a ferromagnet the up spin electron (with the spin 
orientation along the exchange field) decreases its energy by h , while the down spin electron increases its energy by 
the same value. To compensate this energy variation, the up spin electron increases its kinetic energy, while the down 
spin electron decreases its. In the result the Cooper pair acquires a center of mass momentum 2Skp — Ihjvp , which 
implies the modulation of the order parameter with the period -nvp /h. The direction of the modulation wave vector 
must be perpendicular to the interface, because only this orientation is compatible with the uniform order parameter 
in the superconductor. 

To get some idea about the peculiarity of the proximity effect in S/F structures, we may start also from the 
description based on the generalized Ginzburg-Landau functional Eq. (8). Such approach is adequate for a small 
wave- vector modulation case, i. e. in the vicinity of the (H*,T*) point of the (H,T) phase diagram, otherwise the 
microscopical theory must be used. This situation corresponds to a very weak ferromagnet with an extremely small 
exchange field h s=s fisH* = 1.05T C , which is non realistic as usually h >> T c . However, we will discuss this case 
to get a preliminary understanding of the phenomenon. We address the question of the proximity effect for a weak 
ferromagnet described by the generalized Ginzburg-Landau functional Eq. (8). More precisely, we consider the decay 
of the order parameter in the normal phase, i. e. at T > T C i assuming that our system is in contact with another 
superconductor with a higher critical temperature, and the x axis is choosen perpendicular to the interface (see Fig. 
3). 

The induced superconductivity is weak and to deal with it, we may use the linearized equation for the order 
parameter (9), which is written for our geometry as 

fl *-^ + k?= - (12) 

The solutions of this equation in the normal phase are of the type = tp exp(kx), with a complex wave-vector 
k = fci + ift2, and 



*l= { -£ A l + ^^—l , (13) 




*l = ^Z l + \/l + 77^M- (14) 



277 V V T ci - T, 



cu 



If we choose the gauge with the real order parameter in the superconductor, then the solution for the decaying order 
parameter in the ferromagnet is also real 

ip{x) = ip 1 exp(— k\x) cos(fc 2 x), (15) 

where the choice of the root for k is the condition k\ > 0. So the decay of the order parameter is accompanied by 
its oscillation (Fig. 3b), which is the characteristic feature of the proximity effect in the considered system. When 
we approach the critical temperature T ci the decaying wave- vector vanishes, k\ — > 0, while the oscillating wave- vector 



&2 goes to the FFLO wave- vector, fc 2 — ► y so a FFLO phase emerges. Let us compare this behavior with the 
standard proximity effect (Deutscher and De Gennes, 1969) described by the linearized Ginzburg-Landau equation 
for the order parameter 

«^-70 = O, (16) 

with 7 > 0. In such case T c simply coincides with T cu , and the decaying solution is ip = ip exp(— a;/£(T)), where the 
coherence length £(T) = y/j/a (Fig. 3a). This simple analysis brings in evidence the appearance of the oscillations 
of the order parameter in the presence of an exchange field. This is a fundamental difference between the proximity 
effect in S/F and S/N systems, and it is at the origin of many peculiar characteristics of S/F heterostructures. 
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In real fcrromagnets, the exchange field is very large compared with superconducting temperature and energy 
scales, so the gradients of the superconducting order parameter variations are large too, and can not be treated in 
the framework of the generalized Ginzburg-Landau functional. To describe the relevant experimental situation we 
need to use a microscopical approach. The most convenient scheme to do this (see Appendix A and B) is the use 
of the Boboliubov-de Gennes equations or the Green's functions in the framework of the quasiclassical Eilenberger 
(Eilcnberger, 1968) or Usadel (Usadel, 1970) equations. 

If the electron scattering mean free path I is small (which is usually the case in S/F systems), the most natural 
approach is to use the Usadel equations for the Green's functions averaged over the Fermi surface (Appendix). 
Linearized over the pair potential A(x), the Usadel equation for the anomalous function F(x,co) depending only on 
one coordinate x is 

D d 2 \ 

M + ih ■ sgn(u) - Yg^j F{x,u) = A{x), (17) 

where w = (2n + 1) ttT are the Matsubara frequencies, and D = ^vfI is the diffusion coefficient. In the F region, we 
may neglect the Matsubara frequencies compared to the large exchange field (h » T c ), and the pairing potential A 
is absent (we assume that the BCS coupling constant A is zero there). This results in a very simple form of the Usadel 
equation for the anomalous function Ff in the ferromagnet 

D f d 2 F f 

ihsgn(w)F f -^-^f=0, (18) 

where Df is the diffusion coefficient in the ferromagnet. For the geometry in Fig. 3 and lo > 0, the decaying solution 
for Ff is 

F f (x,w > 0) = Aexp ^-li-i^ , (19) 

where £/ = \J is the characteristic length of the superconducting correlations decay (with oscillations) in F- 
layer (see Table I). Due to the condition h >> T c , this length is much smaller than the superconducting coherence 
length £ s = yj 2x7^' i- e - £/ <<; £«■ ^hc cons tant A is determined by the boundary conditions at the S/F interface. 
For example, in the case of a low resistivity of a ferromagnet, at first approximation the anomalous function in a 
superconductor F s is independent on coordinate and practically the same as in the absence of the ferromagnet, i.e. 
F s = A/VA 2 + lo 2 . If, in addition, the interface is transparent then the continuity of the function F at the F/S 
boundary gives A = A/VA 2 + lo 2 . For u> < 0, we simply have Ff (x,ui < 0) = Fj (x,lo > 0). In a ferromagnet, the 
role of the Cooper pair wave function is played by than decays as 

*~]Tf(x,c)~ Aexp(-f )cos(f ). (20) 

We retrieve the damping oscillatory behavior of the order parameter Eq. (15), Fig. 3b. The important conclusion we 
obtain from the microscopic approach is that in the dirty limit the scale for the oscillation and decay of the Cooper 
pair wave function in a ferromagnet is the same. 

In the case of a clean ferromagnet the damped oscillatory behavior of the Cooper pair wave function remains, 
though at zero temperature the damping is non-exponential and much weaker (<~ ^) . Indeed, the decaying solution 
of the Eilenberger equation in the clean limit (see Appendix B) is 

f(x,0,w)~exp[—± '-), 21 

where 9 is the angle between x-axis and Fermi velocity in a ferromagnet, and VFf is its modulus. After averaging 
over the angle and summation over the Matsubara frequencies u we obtain 



7T 

* ~ ^ / f(x,6,uj)sm9d9 ~ -exp( — — )sin(-^-). 

J x ?i/ 6/ 



(22) 
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Here the decaying length £ if — , and the oscillating length £2/ = ^nf ( see Table I) . At low temperature £1 / — ► 
and the Cooper pair wave function decays very slowly ~ i sin(^). An important difference with the proximity effect 
for the normal metal is the presence of the short-ranged oscillations of the order parameter with the temperature 
independent period 2ir{;2 / ■ In contrast with the dirty limit in a clean ferromagnet the characteristic lengths of the 
superconducting correlations' decay and oscillations arc not the same. Halterman and Vails (2001) performed the 
studies of the ferromagnet-superconductor interfaces on the basis of the self-consistent numerical solution of the 
microscopical Bogoliubov-de Gennes equations. They clearly observed the damped oscillatory behavior of the Cooper 
pair wave function of the type * ~ i sin(^). 

We may conclude that at low temperatures the proximity effect in clean ferromagnet metals is long-ranged. On 
the other hand, in the dirty limit the use of the Usadel equations gives the exponential decay of \&. This is due to 
the fact that the Usadel equations are obtained by averaging over the impurities configurations. Zyuzin et al. (2003) 
pointed out that at distances x » £/ the anomalous Green's function F (as well as the Cooper pair wave function) 
has a random sample-specific sign, while the modulus does not decay exponentially. This circumstance leads to the 
survival of the proximity effect in the dirty ferromagnet at distances x » £/. The use of the Usadel equations at 
such distances may be misleading. However, from the practical point of view the range of interest is x < 5£/, because 
at larger distances it is difficult to observe the oscillating phenomena on experiment. In this range the use of the 
Usadel equation is adequate. 

The characteristic length of the induced superconductivity variation in a ferromagnet is small compared with a 
superconducting length, and it implies the use of the microscopic theory of the superconductivity to describe the 
proximity effect in S/F structures. In this context, the calculations of the free energy of S/F structures in the 
framework of the standard Ginzburg-Landau functional (Ryazanov et al, 2001a; Ryazanov et al, 2001b) can not be 
justified. Indeed, the possibility to neglect the higher gradient terms in the Ginzburg-Landau functional implies that 
the length scale of the variation of the order parameter must be larger than the correlation length. In the ferromagnet 

the correlation length is £/ = \J in the dirty limit and £° = -j^- in the clean limit. We see that they coincide with 
the characteristic lengths of the order parameter variation in a ferromagnet. Therefore the higher gradient terms in 
the Ginzburg-Landau functional will be of the same order of magnitude as the term with the first derivative. 



C. Density of states oscillations 

Superconductivity creates a gap in the electronic density of states (DOS) near the Fermi energy Ef, i- e. the DOS 
is zero for a energy E in the interval Ef — A < E < Ef + A. So, it is natural, that the induced superconductivity in 
S/N structures decreases DOS at Ep near the interface. Detailed experimental studies of this phenomenon have been 
performed by Moussy et al. (2001). Damped oscillatory dependence of the Cooper pair wave function in ferromagnet 
hints that a similar damped oscillatory behavior may be expected for the variation of the DOS due to the proximity 
effect. Indeed, the DOS N{e), where e = E — Ep is the energy calculated from the Fermi energy, is directly related 
to the normal Green function in the ferromagnet Gf(x, uj) (Abrikosov et al, 1975) 

Nf (e) = N(0) Re G f (x, lu -> is) , (23) 

where N(0) is the DOS of the ferromagnetic metal. In a dirty limit taking into account the relation between the 
normal and anomalous Green functions + Fj = 1 (Usadel, 1970), and using for Ff = ^ A t +U] 2 CX P ( — 1 £j~ x ) > wc 
directly obtain the DOS at the Fermi energy (e = 0) in a ferromagnet (Buzdin, 2000) at the distance x >> £f 

N f (e = 0) « N(0) (l - Iexp(-^)cos(^)) . (24) 

This simple calculation implies A << T c . An interesting conclusion is that at certain distances the DOS at the Fermi 
energy may be higher than in the absence of superconductor. This contrasts with the proximity effect in the S/N 
systems. Such behavior has been observed experimentally by Kontos et al. (2001) in the measurements of the DOS 
by planar-tunneling spectroscopy in Al/Al 2 03/PdNi/Nb junctions, see Fig. 4. 

For the PdNi layer thickness 50 A we arc at the distance when the term cos(|j) in Eq. (24) is positive and wc have 

the normal decrease of the DOS inside the gap due to the proximity effect. However, for PdNi layer thickness 75 A 
the cos(|j) term changes its sign and the DOS becomes a little bit larger than its value in the normal effect. Such 

inversion of the DOS permits us to roughly estimate £/ for the PdNi alloy used by Kontos et al. (2001) as 60 A. 
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At the moment, there exist only one experimental work on the DOS in S/F systems, while several theoretical papers 
treat this subject more in details. In a series of papers Halterman and Vails (2001, 2002, 2003) performed extensive 
theoretical studies of the local DOS behavior in S/F systems in a clean limit in the framework of the self-consistent 
Bogoliubov-De Gennes approach. They calculated the DOS spectra on both S and F sides and took into account the 
Fermi wave vectors mismatch, interfacial barrier and sample size. 

Fazio and Lucheroni (1999) performed numerical self-consistent calculations of the local DOS in S/F system in the 
framework of the Usadel equation. The influence of the impurity scattering on the DOS oscillations has been studied 
by Baladie and Buzdin, (2001) and Bergeret et al. (2002). An interesting conclusion is that the oscillations disappear 
in the clean limit. In this context it is quite understandable, that the calculations of the DOS oscillations made in 
the ballistic regime for the ferromagnetic film on the top of the superconductor (Zareyan et al., 2001, Zareyan et al., 
2002 ) depend essentially on the boundary conditions at the ferromagnet-vacuum interface. Sun et al. (2002) used the 
quasiclassical version of the Bogoliubov-De Gennes equations for the numerical calculations of the DOS in the S/F 
system with semi-infinite ferromagnet. They obtained in the clean limit the oscillations of the DOS and presented a 
quantitative fit of the experimental data of Kontos et al. (2001). Astonishingly, in the another quasiclassical approach 
on the basis of Eilenberger equations the oscillations of DOS are absent in the case of an infinite electron mean free 
path (Baladie and Buzdin, 2001 and Bergeret et al., 2002). 

DOS oscillations in ferromagncts hint on the similar oscillatory behavior of the local magnetic moment of the 
electrons. The corresponding magnetic moment induced by the proximity effect may be written as 

SM = i(i B N{0)iTTj2 {Gf(x,u,h) - G f (x,cj,-h)) . (25) 
Assuming the low resistivity of a ferromagnet in the dirty limit at temperature near T c , the magnetic moment is 



A 2 9t 9t 

SM = - MB 7V(0)7T— exp(--) sin(-). (26) 
21 c 4/ Kf 



Note that the total electron's magnetic moment in a ferromagnet being 



M = SM + f i B N(0)h. (27) 

Similarly to the DOS the local magnetic moment oscillates, and curiously in some regions it may be higher than in 
the absence of superconductivity. Proximity effect also induces the local magnetic moment in a superconductor near 
the S/F interface at the distance of the order of superconducting coherence length £ s . 

The proximity induced magnetism was studied on the basis of the Usadel equations by Bergeret al. (2004a, 2004b) 
and Krivoruchko and Koshina (2002). Numerical calculations of Krivoruchko and Koshina (2002) revealed the damped 
oscillatory behavior of the local magnetic moment in a superconductor at the scale of £ s with positive magnetization at 
the interface. On the other hand Bergeret al. (2004a) argued that the induced magnetic moment in a superconductor 
must be negative. This is related to the Cooper pairs located in space in such a way that one electron of the pair is in 
superconductor, while the other is in the ferromagnet. The direction along the magnetic moment in the ferromagnet 
is preferable for the electron of the pair located there and this makes the spin of the other electron of the pair (located 
in superconductor) to be antiparallel. 

The microscopic calculations of the local magnetic moment in the pure limit in the framework of Bogoliubov-de 
Gennes equations (Halterman and Vails, 2004) also revealed the damped oscillatory behavior of the local magnetic 
moment but at the atomic length scale. Probably in the quasiclassical approach the oscillations of the local magnetic 
moments disappear in the clean limit, similarly to the case of DOS oscillations. The magnitude of the proximity 
induced magnetic moment is very small, and at present time there are no manifestations of this phenomena on 
experiment. 



D. Andreev reflection at the S/F interface 



The spin effects play an important role in the Andreev reflection at the S/F interface. Indeed, an incident spin up 
electron in ferromagnet is reflected by the interface as a spin down hole, and in the result a Cooper pair of electrons 
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with opposite spins appears in a superconductor. Therefore the both spin up and spin down bands of electrons in 
ferromagnet are involved in this process. De Jong and Bcenakkcr (1995) were the first to demonstrate the major 
influence of spin polarization in ferromagnet on the subgap conductance of the S/F interface. Indeed, in the fully 
spin-polarized metal all carriers have the same spin and Andreev reflection is totally suppressed. In general, with the 
increase of the spin polarization the subgap conductance drops from the double of the normal state conductance to a 
small value for the highly polarized metals. Following de Jong and Beenakker (1995) let us consider a simple intuitive 
picture of the conductance through a ballistic S/F point contact. Using the language of the scattering channels 
(subbands which cross the Fermi level), the conductance at T = of a ferromagnet-normal metal contact is given by 
the Landaucr formula 



Gfn = ^N. (28) 

The total number of scattering channels N is the sum of the spin up and spin down N± channels N = iV-j- +N± , 
and the spin polarization implies that iVj > N± . In the case of the contact of the superconductor with the non- 
polarized metal all electrons are reflected as the holes, which doubles the number of scattering channels and the 
conductance itself. For the spin-polarized metal where iVj > N± , all the spin down electrons will be reflected as the 
spin up holes. However, only the part N± /TVj < 1 of the spin up electrons can be Andreev reflected. The subgap 
conductance of the S/F contact is then 

e 2 / N, \ e 2 

G F s = j(2N l +2iV T ^)=4-iV i . (29) 

Comparing this expression with Eq. (28) we see that Gfs/Gfn — 4A^ /(N± + ) < 2 and Gps = for the 
full-polarized ferromagnet with N± = 0. If the spin polarization is defined as P = (iVj — iVj_ ) / (ATj_ + ATj ), then the 
suppression of the normalized zero-bias conductance gives the direct access to the value of P : 



_£i = 2(l-P). (30) 

UFN 

The subsequent experimental measurements of the spin polarization with Andreev reflection (Upadhyay et al., 1998; 
and Soulen et al, 1998 J fully confirmed the efficiency of this method to probe the ferromagnets. The Andreev point 
contact spectroscopy permits to measure the spin polarization in a much wider range of materials (Zutic, Fabian and 
Das Sarma, 2004) comparing with the spin-polarized electron tunneling (Meservey and Tedrow, 1994). 

However, the interpretation of the Andreev reflection data on the conductance of the S/F interfaces and the 
comparison of the spin polarization with the one obtained from the tunneling data, may be complicated by the band 
structure effects (Mazin, 1999). Zutic and Vails (1999, 2000), Zutic and Das Sarma (1999) generalized the results of 
the theoretical analysis of Blonder, Tinkham and Klapwijk (1982) to the case of the S/F interface. An interesting 
striking result is that in the absence of the potential barrier at the S/F interface, the spin polarization could increase 
the subgap conductance. The condition of perfect transparency of the interface is Vf^Vfi = w 2 , where Vfi and Vf\, 
are the Fermi velocities for two spin polarizations in ferromagnet, and v s is the Fermi velocity in superconductor. 
Vodopyanov and Tagirov (2003a) proposed a quasiclassical theory of Andreev reflection in F/S nanocontacts and 
analyzed the spin polarization calculated from the conductance and tunneling measurements. 

Note that a rather high spin polarization has been measured in CrC>2 films P = 90% and in Lao.7Sro.3Mn03 films 
P = 78% (Soulen et at, 1998). The spin-polarized tunneling data for these systems is lacking. 

Another interesting effect related with the crossed Andreev reflection has been predicted by Dcutshcr and Feinberg 
(2000) (see also Deutsher, 2004 and Yamashita, Takahashi and Maekawa, 2003). The electric current between two 
ferromagnetic leads attached to the superconductor strongly depends on the relative orientation of the magnetization 
in these leads. If we assume that the leads are fully polarized, then the electron coming from one lead cannot 
experience the Andreev reflection in the same lead. However, this reflection is possible in the second lead, provided its 
polarization is opposite, and the distance between the leads is smaller than the superconducting coherence length. The 
resistance between the leads will be high for the parallel orientation of the magnetizations and low for the antiparallcl 
orientation. 
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IV. OSCILLATORY SUPERCONDUCTING TRANSITION TEMPERATURE IN S/F MULTILAYERS AND BILAYERS 



A. First experimental evidences of the anomalous proximity effect in S/F systems 

The damped oscillatory behavior of the superconducting order parameter in ferromagnets may produce the com- 
mensurability effects between the period of the order parameter oscillation (which is of the order of £/) and the 
thickness of a F layer. This results in the striking non-monotonous superconducting transition temperature depen- 
dence on the F layer thickness in S/F multilayers and bilayers. Indeed, for a F layer thickness smaller than , the 
pair wave function in the F layer changes a little and the superconducting order parameter in the adjacent S layers 
must be the same. The phase difference between the superconducting order parameters in the S layers is absent and 
we call this state the "0"-phase. On the other hand, if the F layer thickness becomes of the order of £/, the pair 
wave function may go trough zero at the center of F layer providing the state with the opposite sign (or tt shift of the 
phase) of the superconducting order parameter in the adjacent S layers, which we call the 'V-phase. The increase of 
the thickness of the F layers may provoke the subsequent transitions from "0"- to "7r" -phases, what superpose on the 
commcnsurability effect and result in a very special dependence of the critical temperature on the F layer thickness. 
For the S/F bilayers, the transitions between "0" and "7r" -phases are impossible; the commensurability effect between 
£f and F layer thickness nevertheless leads to the non- monotonous dependence of T c on the F layer thickness. 

The predicted oscillatory type dependence of the critical temperature (Buzdin and Kuprianov, 1990; Radovic et 
al, 1991) was subsequently observed experimentally in Nb/Gd (Jiang et al, 1995), Nb/CuMn (Mercaldo et al, 1996) 
and Nb/Co and V/Co (Obi et al, 1999) multilayers, as well as in bilayers Nb/Ni (Sidorenko et al, 2003), trilayers 
Fc/V/Fc (Garifullin et al, 2002), Fe/Nb/Fe (Miihge et al, 1996), Nb/[Fc/Cu] layers (Vclcz et al, 1999) and Fc/Pb/Fc 
(Lazar et al, 2000). 

The strong pair-breaking influence of the ferromagnet and the nanoscopic range of the oscillation period complicate 
the observation of this effect. Advances in thin film processing techniques were crucial for the study of this subtle 
phenomenon. The first indications on the non-monotonous variation of T c versus the thickness of the F layer was 
obtained by Wong et al (1986) for V/Fe superlattices. However, in the subsequent experiments of Koorevaar et 
al (1994), no oscillatory behavior of T c was found, while the recent studies by Garifullin et al (2002) of the 
superconducting properties of Fe/V/Fe trilayers even revealed the re-entrant T c behavior as a function of the F layer 
thickness. Bourgeois and Dynes (2002) studied amorphous Pb/Ni bilayer quench-condensed films and observed only 
monotonic depairing effect with the increase of the Ni layer thickness. In the work of Sidorenko et al. (2003), the 
comparative analysis of different techniques of the sample preparation was made and the conclusion is, that the 
molecular beam epitaxy (MBE) grown samples do not reveal T c oscillations, whereas magnetron sputtered samples 
do. This difference is attributed to the appearance of magnetically "dead" interdiffused layer at the S/F interface 
which plays an important role for the MBE grown samples. The transition metal ferromagnets, such as Fc, have a 
strongly itinerant character of the magnetic moment which is very sensitive to the local coordination. In thin Fe 
layers, the magnetism may be strongly decreased and even vanished. Probably the best choice is to use the rare-earth 
ferromagnetic metal with localized magnetic moments. This has been done by Jiang et al (1995) who prepared the 
magnetron sputtered Nb/Gd multilayers, which clearly revealed the T c oscillations, Fig. 5. 

The curves show a pronounced non-monotoneous dependence of T c on the Gd layer thickness. The increase of T c 
implies the transition from the "0" -phase to the "7r" -phase. Note that the previous experiments on the MBE grown 
Nb/Gd samples (Strunk et al, 1994) only revealed the step-like decrease of T c with increasing Gd layer thickness. 
The comprehensive analysis of different problems related to the samples quality was made by Chien and Reich (1999). 
Aarts et al (1997), studied in detail the proximity effect in the system consisting of the superconducting V and 
ferromagnetic Vi- x Fe x alloys and demonstrated the important role of the interface transparency for the understanding 
of the pair-breaking mechanism. 

B. Theoretical description of the S/F multilayers 

To provide the theoretical description of the non-monotoneous dependence of T c , we consider the S/F multilayered 
system with a thickness of the F layer 2df and the S layer 2d s , see Fig. 6. 

The a;-axis is chosen perpendicular to the layers with x — at the center of the S layer. The "0" -phase case 
corresponds to the same superconducting order parameter sign in all S layers (Fig. 6a) while in the "7r" -phase the 
sign of the superconducting order parameter in adjacent S layers is opposite (Fig. 6b). In the case of a S/F bilayer, 
the anomalous Green function F{x) has zero derivative at the boundary with vacuum, see Eq. (32) below. It is 
just the case for the function F(x) in the "0"-phase at the centers of the S and F layers. So the superconducting 
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characteristics of a S/F bilayer with thicknesses d s and df of the S and F layers respectively are equivalent to that of 
the S/F multilayer with double layer thicknesses (2d s and 2df). 

The approach based on the quasiclassical Eilenberger (1968) or Usadel (1970) equations is very convenient to deal 
with S/F systems (see Appendix B). In fact, it is much simpler than the complete microscopical theory, it does not 
need the detailed knowledge of all the characteristics of the S and F metals, and is applicable for scales larger than 
the atomic one. Then, it must work for thicknesses of the layers in the range 20 — 200 A , which is of primary interest 
for S/F systems. 

In the dirty limit, if the electron elastic scattering time r = l/vj is small, more precisely T c t <C 1 and hr <C 1, 
the use of the Usadel equations is justified. The second condition, however is much more restrictive due to a large 
value of the exchange field (h ^> T c ). The Usadel equations deal only with the Green's functions G(x, io) and F(x, u>) 
averaged over the Fermi surface. Moreover, to calculate the critical temperature of the second-order superconducting 
transition in S/F systems, it is enough to deal with the limit of the small superconducting order parameter (A — > 0) 
in the Usadel equations. This linearization permits to put G — sgn(uj) and in the form linearized over A, the Usadel 
equation for the anomalous function F 8 in the S region is written as 

»\ = ( 31 ) 



2 Ox 2 , 

where D s is the diffusion coefficient in the S layer. In the F region, the exchange field is present while the pairing 
potential A is absent, and the corresponding Usadel equation for the anomalous function Ff is just the Eq. (18). 

The equations for F s and Ff must be supplemented by the boundary conditions. At the superconductor- vacuum 
interface, the boundary condition is simply a zero derivative of the anomalous Green function, which implies the 
absence of the superconducting current through the interface. The general boundary conditions for the Usadel 
equations at the superconductor-normal metal interface have been derived by Kupriyanov and Lukichev (1988) and 
near the critical temperature they read 



V dx J x=0 <j s V dx J x=Q 



F s (0) = F f (0) - UlB ( ) , (32) 



x=0 



where 07 (a s ) is the conductivity of the F-layer (S-layer above T c ) . The parameter 75 characterizes the interface 
transparency T = pr~ and is related to the S/F boundary resistance per unit area Ri, via the following simple 

relationship jb — R ^ af (Kupriyanov and Lukichev, 1988). By analogy with the superconducting coherence length 

£ s = \J 2^t ' we introduce the normal metal coherence length £„ = \J . The presented form of the boundary 
conditions corresponds to the S/F interface x = and the positive direction of the x axis chosen along the outer normal 
to the S surface (i. e. the x axis is directed from the S to the F metal). It is worth notify that the boundary conditions 
for the Usadel equations (Kupriyanov and Lukichev, 1988) have been obtained for superconductor/normal metal 
interfaces, and their applicability for S/F interfaces is justified, provided that the exchange field in the ferromagnet 
is much smaller than the Fermi energy, i. e. h << Ep. For a ferromagnet with localized moments, such as Gd, 
this condition is always fulfilled, while it becomes more stringent for transition metals and violated for half-metals. 
Recently Vodopyanov and Tagirov (2003b) obtained the boundary conditions for Eilenberger equations in the case of 
a strong ferromagnet. They used them to study the critical temperature of a S/F bilayer when ferromagnet is in the 
clean limit. Nevertheless the important question about the boundary conditions for Usadel equations at the interface 
superconductor/strong ferromagnet is still open. 

Provided the solutions of Usadel equations in the F and S layers are known, the critical temperature T* may be 
found from the self-consistency equation for the pair potential A(x) in a superconducting layer 

A(x)=irT:\Y / Fs(x,u;), (33) 

where A is BCS coupling constant in S layer (while in F layer it is supposed to be equal to zero). This equation is 
more convenient to write in the following form 

A^ln^+Tr^W^-i^a^)) =0, (34) 
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where T c is the bare transition temperature of the superconducting layer in the absence of the proximity effect. 

The Usadel equations provide a good basis for the complete numerical solution of the problem of the transition 
temperature of S/F superlattices. Firstly such a solution has been obtained for a S/F system with no interface barrier 
by Radovic et al. (1988, 1991), using the Fourier transform method, and this case was treated analytically by Buzdin 
and Kuprianov (1990) and Buzdin et al. (1992). The role of the S/F interface transparency has been elucidated by 
Proshin and Khusainov (1997), (for more references see also the review by Izyumov et al. 2002) and Tagirov (1998). 
Recently Fominov et al. (2002), performed a detailed analysis of the non-monotonous critical temperature dependence 
of S/F bilayers for arbitrary interface transparency and compared the results of different approximations with exact 
numerical calculations. 

Below we illustrate the appearance of the non-monotonous superconducting transition temperature dependence for 
the case of a thin S-layer, which has a simple analytical solution. More precisely, we consider the case d s <C £ s , which 
implies that the variations of the superconducting order parameter and anomalous Green's function in the S layer are 
small. We may write the following expansion up to the x 2 order term for F s in the S layer centered at x = : 

F s (x,w) = F (l-^x 2 y (35) 

where Fq is the value of the anomalous Green's functions at the center of the S-layer, and the linear over x term is 
absent due to the symmetry of the problem in both "0"- and 'V -phases (see Fig. 4). Putting this form of F s into 
the Usadel equation (31), we readily find 

= — ^rr, (36) 

LU + T S 

where we have introduced the complex pair-breaking parameter r" 1 = ^-/3 U and in the first approximation over 
d s /£ s <C 1, the pair potential A may be considered as spatially independent. The pair-breaking parameter r" 1 , is 
directly related to the logarithmic derivative of F s at x — d s 

The boundary conditions Eq. (32) permit us to calculate the parameter t~ x , provided the anomalous Green function 
in the F layer is known: 

!_ DsVf F' f (d s )/F f (d s ) 

2d s a s l-Z nlB F' } (d s )/F f {d s y 1 ' 



C. "0"- and "7r"-phases 

The solution of the Usadel equation (18) in the F layer is straightforward but different for "0"- and "7r"-phases. 
Let us start first with a "0" -phase. In such a case (see Fig. 6a), we must take as a solution for Ff(x) at lo > in the 
interval d s < x < d s + 2df the function symmetrical relative to the plane x = d s + df, i. e. 



Ff(x,u> > 0) = A cosh 



i + 1 



(x-d s - d f ) 



(39) 



Therefore the pair-breaking parameter r^ 1 for "0" -phase at lo > is 



D s aj_ £+1 tanh (jjrj/ 

2d s a s i f 1+ i±l^ Btanh (^i±l df y 



(40) 



and does not depend on the Matsubara frequencies lo. For a negative lo we simply have r s q (lo < 0) = (r s q (lo > 0))*. 
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Now, let us address the case of the "tt" -phase. The only difference is that in such case we must choose the 
asymmetrical solution for Ff(x) 



F f (x,u> > 0) = Bsinh 



i + 1 




(x - d s - df) 



and the corresponding pair-breaking parameter t s * is given by the expression 



7--i(w>0) = (r-i(a;<0))* = 



(41) 



(42) 



D s crj i + 1 



coth 



24 t s ^ j + i±!£„ 7B coth ( ^fd/ 



We see that in all cases the pair-breaking parameter t" 1 is complex and depends on the sign of the Matsubara 
frequency only but not on its value. As a result, with the help of the self-consistency equation (34), we obtain the 
following expression for the critical temperature T* of the S/F multilayer 



ln^ 



*'5 



1 



2ttT*t s 



(43) 



where \& is the Digamma function, and the pair-breaking parameter r s _1 is given by Eqs. (40) and (42) for the "0"- and 
"7r" -phases respectively. This type of expression for T* reminds the corresponding formula for the critical temperature 
of a superconductor with magnetic impurities (Abrikosov and Gor'kov, 1960), though the "magnetic scattering time" 

'" << 1), the formula for the critical 



t s is complex in our system. If the critical temperature variation is small ( — ^ 
temperature shift Eq. (43) may be simplified 



iP—k"^- 



(44) 



D. Oscillating critical temperature 

To illustrate the oscillatory behavior of the critical temperature, we consider the case of a transparent S /F interface 
7s = 0. The critical temperatures T*° and T* 71 for the "0"- and "7r" -phases respectively, are 



T r - T, 



*o 



/ sinh(2y) - sin(2y) \ 

^ cosh(22/) + cos(2y) ) ' 1 °> 



T c AT c t Vcosh(2y) + cos(2y) J ' 

T c - T™ _ tt f sinh(2y) + sin(2y) \ 



4T c r Vcosh(2y) - cos(2y) / ' 
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where t 1 = 2 d\ } V~ ana - 2j/ = ^df/£f is the dimensionless thickness of the F layer. The critical temperature variation 
versus the F layer thickness is presented in Fig. 7. 

We see that for the small F layer thicknesses, the "0" -phase has a higher transition temperature. The first crossing 
of the curves T*°(y) and T c * 7r (y) occurs at 2y c « 2.36 and in the interval of thickness 2.36£/ < 2df < 5.5£/ , the 'V'- 
phase has a higher critical temperature. The oscillations of the critical temperature rapidly decay with the increase 
of y, and it is not realistic to observe on experiment more than two periods of oscillations. 

In the general case, the F-layer thickness dependence of the critical temperature Eq. (43) may be written for 
"0" -phase in the following form convenient for numerical calculations 

T*° /1\ 

'"ir=*U)- (47) 

,.- . ,' . -IT,. 
-Re* 



T*°f 7 + ^ coth [(1 + i)y] 



where the dimensionless parameter f 1 = 1/ (4ttT c to) and 7 = 7s (&»/£/) ■ The corresponding formula for the critical 
temperature for the "tt" -phase is simply obtained from Eq. (47) by the substitution coth — > tanh. 
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In Fig. 8, we present the examples of calculations of the thickness dependence of the critical temperature for S/F 
multilayers for different interface transparencies. 

The oscillations of the critical temperature are most pronounced for transparent interface 7 = 0, and they rapidly 
decrease with the increase of the boundary barrier (at 7 > 2 the oscillations are hardly observable). Note that, for 
certain values of the parameters fo and 7, the T*°(df) dependence may show the infinite derivative, which indicates 
the change of the order of the superconducting transition from second-order to the first-order one. This question was 
studied in detail by Tollis (2004) . The increase of the boundary barrier not only decreases the amplitude of the critical 
temperature oscillations, but also it decreases the critical thickness of F layer y c , corresponding to the "0"- 'V-phase 
transition. The limit 7 = 75 (£«/£/) >> 1 i s rat her special one. In such case the S/F interface barrier becomes a 
tunnel barrier, and the critical thickness y c may be much smaller than 1. Indeed, if the critical temperature variation 
is small, (more precisely if jt » 1), the condition Re(r~J) = Re(r~ 7 J) is realized at 



and the mechanism of the "0"- "it" -phase transition is now related to the peculiarity of tunneling through the F 
layer. This is very different from the case of low interface transparency, when the transition occurs due to the spatial 
oscillations of the anomalous Green's function. It must be very difficult to observe the low transparency regime of 
the "0"- 'V transition with the help of the critical temperature measurements due to the fact that at 7 >> 1 the 
oscillations of T*(df) become very small. On the other hand, the measurements of the critical current in S/F/S 
Josephson junctions may be the adequate technique to reveal the "0"- "7r" transition in this regime (see next Section). 

It is interesting to note that for small thicknesses of F layer (df < £/) the critical temperature decreases with the 
increase of the interface barrier (provided the condition j(df/£f) < 1 is fulfilled) - see Fig. 8. Such a counterintiutive 
behavior may be explained in the following way. The low penetration of the barrier prevents the quick return of the 
Cooper pair from thin F layer. Therefore, the Cooper pair stays for a relatively long time in the F layer before going 
back to the S layer. In the results, the pair-breaking role of the exchange field in the F layer occurs to be strongly 
enhanced. 

The cases of S/F bilayers or F/S/F threelayers with parallel magnetization are equivalent to the "0" -phase case 
for the multilayers (with double F layers thickness) and the corresponding T* (df) dependence reveals a rather weak 
non-monotonous behavior in the case of finite transparency of the S/F interface (see Fig. 8). The comparison of 
the experimental data of Ryazanov et al. (2003) for the critical temperature of the bilayer Nb/Cuo. 43^0.57 vs the 
thickness of the ferromagnetic layer with the theoretical fit (Fominov et al., 2002) is presented in Fig. 9. 

Now let us address a question, if it is possible to have a transition into a state with the phase difference another 
than and it ? For example the state with the intermediate phase difference < ipo < n may be expected at F 
layer thicknesses near dj. The numerical calculations of Radovic et al. (1991) indeed revealed the presence of the 
intermediate phase. However, the relative width of the region of its existence near dfj was very small - around several 
percents only. On the other hand, the analytical calculations show that for the thin S layer case the states without 
current (corresponding to the highest T*) are possible only for the phase difference or n . Also, in the S/F/S 
junctions the transitions between "0"- and " 7r" -states are discontinuous - see discussion in the next section. Probably 
the narrow region of the "ipo"- phase existence obtained by the numerical calculations (Radovic et al. , 1991) is 
simply related with its accuracy ~ 1%, and the width of this region may decrease with the increase of the accuracy. 
Nevertheless there is another mechanism of the realization of the "<po"- phase due to the fluctuations of the thickness 
of F layer. In such case near the critical F layer thickness d C f the regions of "0"- and " 7r" -phases would coexist. If the 
characteristic dimensions of these regions are smaller than the Josephson length in S/F structure, then the average 
phase difference would be different from and 7r (Buzdin and Koshclcv, 2003). 

The quasiclassical Eilenberger and Usadel equations are not adequate for treating the strong ferromagnets with 
h ~ Ep because the period of Green's function oscillations becomes comparable with the interatomic distance. On 
the other hand, the approach based on the Bogoliubov-dc Gcnncs equations in clean limit is universal. Halterman and 
Vails (2003, 2004a) applied it to study the properties of clean S/F multilayers, at low temperature. They obtained the 
excitation spectrum through numerical solution of the self-consistent Bogoliubov-De Gennes equations and discussed 
the influence of the interface barrier and Fermi energy mismatch on the local density of states. Comparing the energy 
of the "0" and "7r" phases Halterman and Vails confirmed the existence of the transitions between them with the 
increase of F layer thickness. It is of interest that the local density of states is quite different in the "0" and "it" 
phases, and its measurements could permit to trace the "0" - "7r" transition. In the more recent work Halterman 
and Vails (2004b) showed that a lot of different order parameter configurations may correspond to the local energy 
minima in S/F heterostructures. 
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The calculations of the energy spectrum in the S/F/S system in "0" and "tt" phases on the basis of Eilenberger 
equations were performed by Dobrosavljevic-Grujic, Zikic and Radovic (2000) for s-wave superconductivity and d- 
wave superconductivity (Zikic et al, 1999). The large peaks in the density of states were attributed to the spin-split 
bound states appearing due to the special case of the Andreev reflection at the ferromagnetic barrier. 

In the previous analysis the spin-orbit and magnetic scattering were ignored. Dcmlcr, Arnold, and Bcasley (1997) 
theoretically studied the influence of the spin-orbit scattering on the properties of S/F systems and demonstrated 
that it is quite harmful for the observation of the oscillatory effects. A similar effect is produced by the magnetic 
scattering which at some extend is always present in S/F systems due to the non-stoichiometry of the F layers (and 
it may be rather large when the magnetic alloy is used as F layer). The calculations of the critical temperature 
of the S/F multilayers in the presence of the magnetic scattering were firstly performed by Tagirov (1998). In the 
framework of the formalism presented in this section it is very easy to take into account the magnetic diffusion with 
the spin-flip scattering time r m - it is enough to substitute the exchange field h in the linearized Usadel equation 
(17) by h — js5n(w)T~ 1 . This renormalization leads to the decrease of the damping length and the increase of the 
oscillation period, which makes the T*(dj) oscillations less pronounced (Tagirov, 1998). 

V. SUPERCONDUCTOR-FERROMAGNET-SUPERCONDUCTOR 'V-JUNCTION 

A. General characteristics of the "n" -junction 

A Josephson junction at equilibrium has usually a zero phase difference tp between two superconductors. The energy 
E of Josephson junction may be written as (see for example De Gennes, 1966a) 

£7=|^(l-cos V ), (49) 

where I c is the Josephson critical current, and the current-phase relation is I s (<p) — ^ ^ = I c simp. At the standard 
situation, the constant I c > 0, and the minimum energy of a Josephson junction is achieved at ip = 0. However, 
in the previous section it has been demonstrated that in the S/F multilayers the transition into the "tt" -phase may 
occur. This means that for the Josephson S/F/S junction (with the same thickness of F layer which corresponds to 
the "7r" -phase in the multilayered system) the equilibrium phase difference would be equal to tt, and it is natural to 
call such a junction the "tt" -junction. For the "tt" -junction, the constant I c in the equation (49) is negative, and 
the transition from "0"- to "7r" -state may be considered as a change of the sign of the critical current, though the 
experimentally measured critical current is always positive and equals to |/ c | . The S/F/S junctions would reveal the 
striking non-monotonous behavior of the critical current as a function of F layer thickness. The vanishing of the 
critical current signals the transition from "0"- to "tt" -state. 

The possibility of the negative Josephson coupling was firstly noted by Kulik (1966), who discussed the spin-flip 
tunneling through an insulator with magnetic impurities. Bulaevskii et al. (1977) put forward the arguments that 
under certain conditions such a spin- flip tunneling could dominate the direct tunneling and lead to the "7r" -junction 
appearance. Up to now there are no experimental evidences of the " tt" -coupling in the Josephson junctions with 
magnetic impurities. On the other hand, Buzdin et al. (1982) showed that in the ballistic S/F/S weak link I c displays 
damped oscillations as a function of the thickness of the F layer and its exchange field. Later, Buzdin and Kuprianov 
(1991) demonstrated that these oscillations remain in the diffusive regime and so, the " tt" -coupling is the inherent 
property of the S/F/S junctions. The characteristic thickness of F layer corresponding to the transition from the "0"- 

to "7r" -phase is = \J ^f, and it is rather small (10 — 50) A in the typical ferromagnets because of the large value 
of the exchange field (h > lOOOif). So, the experimental verification of the "tt" -coupling in S/F/S junction was not 
easy, due to the needed very careful control of the F layer thickness. Finally the first experimental evidence for a "tt"- 
junction was obtained by Ryazanov et al. (2001a) for the Josephson junction with a weakly ferromagnetic interlayer 
of a Cu x Nii_ x alloy. Such choice of F layer permitted to have a ferromagnet with a relatively weak exchange field 
(h ~ 100 — 500 K) and, therefore the relatively large £/ length. 

B. Theory of " tt" -junction 

The complete qualitative analysis of the S/F/S junctions is rather complicated, because the ferromagnetic layer 
may strongly modify superconductivity near the S/F interface. In addition, the boundary transparence and electron 
mean free path, as well as magnetic and spin-orbit scattering, are important parameters affecting the critical current. 
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To introduce the physics of "7r" -coupling, we prefer to concentrate on the rather simple approach based on the 
Usadel equation and consider the S/F/S junction with a F-layer of thickness 2df, (see Fig. 10) and identical S/F 
interfaces. In the case of small conductivity of F layer or small interface transparency cr/£ s /<7 s £/ « max(l,-f B ) we 
may use the "rigid boundary" conditions (Golubov et al, 2004) with F s (-d f ) = Ae^/ 2 /V^ 2 + A 2 and F s (d f ) = 
Ae^/ 2 /V^ 2 + A 2 . 

The solution of Eq. (18) in a fcrromagnet satisfying the corresponding boundary conditions is written as 



F(x) = 



Vto 2 + A 2 

cos (ip/2) cosh (kx) 
(cosh (kdf) + kj B t; n sinh (kdf)) 
i sin (tp/2) sinh (kx) 
(sinh (kdf) + k^ B ^ n cosh (kdf)) 



(50) 



where the complex wave-vector k = y/2 (juJj + isign (uu) h) /Df. This solution describes the F(x) behavior near the 
critical temperature. Note, that in principle, at arbitrary temperature, the boundary conditions are different from the 
Eq. (32), see for example (Golubov et al., 2004). However, in the limit of low S/F interface transparency ("f B » 1), 
when the amplitude of the F function in F- layer is small, we may use the linearized Usadel equation (18) at all 
temperatures. The only modification in the boundary conditions Eq. (32) is that F s must be substituted by F s / \ G S \ 
and 7 B by 7 B / |G S |, where the normal Green function in superconducting electrode G s = cu/y/oj 2 + A 2 . Taking this 
renormalization into account in the explicit form Eq. (50) , we may use it in the formula for the supercurrent 



I s ( V ) = ieN(0)Df*T SJ2[F^F-F-^F 



dx 



(51) 



where F (x, h) — F*(x, —h), S is the area of the cross section of the junction and N(0) is the electron density of state 
for a one spin projection. This expression gives the usual sinusoidal current-phase dependence I s (ip) = I c sin(tp) with 
the critical current 



I c = eSN(0)D f TTT 

A 2 2k/cosh(2kd f ) 



E 



uj 2 tanh (2kd f ) (1 + T 2 fc 2 ) + 2kT u 



(52) 



This expression may be easily generalized to take into account the different interface 

2 . I<- I \n l\2 



where T w = 7 B 6»/I G . 

transparencies j B1 , r ) B2 » 1, it is enough to substitute in Eq. (52) T 2 — ► 7 B1 7 B2 (£ n / \G S \) Z and 2r 
{Ibi + I132) £n/ \G S \. Near T c and in the case of transparent interface j B — > (Buzdin and Kuprianov, 1991) 



I c = eSN(0)D f 



ttA 2 
~2Tc 



Re 



sinh (2kdf) 



(53) 



tin 



cos (2y) sinh (2y) + sin (2y) cosh (2y) 



cosh (Ay) — cos (Ay) 



where 2y — 2df/£f is the dimensionless thickness of the F layer, R n = 2df / (<JfS) is the resistance of the junction 
(a f = 2e 2 N(0)D f is the conductivity of the F layer), and V = f^. 

The dependence I c R n /Vo vs. 2y is presented in Fig. 11. The first vanishing of the critical current signals the 
transition from "0"- to "7r" -state. It occurs at 2y c « 2.36 which is exactly the critical value of F layer thickness 
in S/F multilayer system corresponding to the "0"— "tt" -state transition, i. e. to the condition T*° = T* v in the 
Eqs.( 45). The theoretical description of the S/F/S junctions with arbitrary interface transparencies near the critical 
temperature was proposed by Buzdin and Baladie (2003). 

At low temperature or low S/F barrier the amplitude of the anomalous Green's function Ff(x) is not small and we 
need to use the complete (non-linearized) Usadel equation. In the limit of large thickness of F layer df » £/ and 
7 B = 0, the analytical solution was obtained by Buzdin and Kuprianov (1991), and the critical current is 



I c R n = 64V2^^ (^Pj 2ycxp(-2y) |sin fey + J) 



(54) 
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with the function 



^)=^£ - 2 , (55) 



20 + + w 



where f2 = \J uj 2 + |A| 2 , and J 7 (t^) ~ ifs ^ at T « T c while at low temperature T << T Cl the function (^-^ ~ 
0.071. 

Note that in the clean limit (rh » 1) the thickness dependence of the critical current is very different (Buzdin et 
al, 1982) and near T c it is 



IcRn — 



■ I ihdt \ 



4e ( ihd f ) 

V Up I 



(56) 



i. e. the critical current decreases ^ 1/df and not exponentially like in the dirty limit case. In general, in the clean 
limit the S/F proximity effect is not exponential, but a power low one. 

The expression (56) was obtained on the basis of Eilenberger equations. In the case of a strong ferromagnet 
h < Ef , the period of the oscillations of the Green's functions becomes of the order of the interatomic distance, 
and this approach does not work anymore. Using the technique of the Bogoliubov-de Gennes equations, Cayssol and 
Montambaux (2004) demonstrated that the quasiclassical result (56), where the only relevant parameter for the critical 
current oscillations being hdf/vF, is not applicable for the strong ferromagnets. This is related to the progressive 
suppression of the Andreev reflection channels with the increase of the exchange energy. 

In the framework of the Bogoliubov-de Gennes equations Radovic et al. (2003) studied the general case of the bal- 
listic S/F/S junction for a strong exchange field, arbitrary interfacial transparency and Fermi wave vectors mismatch. 
The characteristic feature of such ballistic junction is the short-period geometrical oscillations of the supercurrent as 
the function of df due to the quasiparticle transmission resonances. In the case of strong ferromagnet, the period of 
"0"— "tt" oscillations becomes comparable with the period of geometrical oscillations, and their interplay provides 
very special I c (df) dependences. Also Radovic et al. (2003) demonstrated that the current-phase relationship may 
strongly deviate from the simple sinusoidal one, and studied how it depends on the junction parameters. While the 
temperature variation of I c is usually a monotonic decay with increasing temperature, near the critical thickness df 
corresponding to "0"— "7r" transition, a nonmonotonic dependence I c on temperature was obtained. Radovic et al. 
(2001) showed that at low temperature the characteristic multimodc anharmonicity of the current-phase relation in 
clean S/F/S junctions implies the coexistence of stable and metastable "0—" and "ir— " states. As a consequence, the 
coexistence of integer and half-integer flixoid configuration of SQUID was predicted. Note that for strong ferromagnets 
the details of the electrons energy bands become important for the description of the properties of S/F/S junction. 

The weak link between d— wave superconductors may also produce the it shift effect (as a review, see for exam- 
ple Van Harlingcn, 1995). The situation of the Joscphson coupling in a ferromagnetic weak link between d— wave 
superconductors was studied in the clean limit theoretically by Radovic et al. (1999). 

It is interesting that in the limit kd << 1 (i.e. df << £f ) the oscillations of the anomalous function in the F 
layer are absent, but as it has been noted previously, for the case of the low transparency of the barrier 7 B >> 1, the 
critical current can nevertheless change its sign. Indeed, in this limit, the expression for the critical current Eq. (52) 
reads 



2 IAI 1 / 1 

I c = eN(0)D f irT S V !— ' 

2rf/ 2 1 M 



lsindfk\: , 2 + | A |2 



(57) 



Usually at experiment, the Curie temperature O of ferromagnet is higher than the superconducting critical temperature 
T c . For RKKY mechanism of ferromagnetic transition ~ h 2 /Ep and so the exchange field h occurs to be much 
larger than the superconducting critical temperature T c . In the case of the itinerant ferromagnetism, the exchange 
field is usually several times higher than the Curie temperature and also the limit h >> T c holds. Taking this into 
account and performing the summation over Matsubara frequencies of the first two terms in the brackets of the Eq. 
(57), we finally obtain 
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_ eN(0)SD f Ae f 
A 



h 



+ (58) 



2ttTA0 ^ lu 4tt d 2 f \ ( A\ 

T7 22 71 , 7^372 -T bf tanh 



7 B Cnrf/ ^ (^ + A2) 3/2 3 V217J- 

We start with the analysis of I c over c?/ dependence in the limit of very large 7 B (more precisely when j B » (r) )■ 
In such case we may neglect the term proportional to 1/7 B in the brackets of Eq. (58), and then we obtain that at 
T — > the transition into the 7r— phase occurs ( I c changes its sign) at 



Indeed the condition df « is satisfied. In the case of very low boundary transparencies, the relevant formula 
obtained in (Buzdin and Baladie, 2003) near the critical temperature in the limit (T c /h) — ► also reveals the 
crossover between "0"— and "7r"— phase. On the other hand, no transition into 7r— phase was obtained in the analysis 
of S/F/S system by Golubov et al. (2002b), which is apparently related to the use of the gradient expansion of the 
anomalous function in ferromagnet when only the first term has been retained. 

It is interesting to note that the critical F-layer thickness d^, when the transition from "0"— to "7r"— phase occurs, 
depends on the temperature. The corresponding temperature dependences are presented in Fig. 12 for different value 
of (T c /h) ratios. We see that dj(T) decreases when the temperature decreases. This is a very general feature and it is 
true also for the subsequent "0" — "7r" transitions occurring at higher F layer thickness. So for some range of F-layer 
thicknesses the transition from "0"— to "tt"— phase is possible when the temperature lowers. 

For the case of moderately large j B , i.e. when 1 « j B « the terms with \P functions in Eq. (58) can be 
neglected, and at T = T c the critical thickness dfj is 

d c f (T = T c ) = ^ (^L) ^ , (60) 



2 



1/3 



while at T — > the critical thickness is somewhat smaller dj (T = 0) = 4f J ■ The critical F layer thickness, 

given by Eq. (60), naturally coincides with the corresponding expression Eq. (48) obtained for S/F multilayers in 
the limit h >> T c . The examples of different non- monotonous I C (T) dependences for low barrier transparency limit 

7 B >> (jj^j are presented in Fig. 13. In fact, in the limit of low barrier transparency and thin F layer, we deal 
with the superconducting electrons tunneling through ferromagnetically ordered atoms. The situation is in some 
sense reminiscent the tunneling through magnetic impurities, considered by Kulik (1966) and Bulacvskii et al. (1977). 
What may be more relevant is the analogy with the mechanism of the 'V— phase realization due to the tunneling 
through a ferromagnetic layer in the atomic S/F multilayer structure, which we consider in the section 7. 

Fogelstrom (2000) considered the ferromagnetic layer as a partially transparent barrier with different transmission 
for two spin projections. In some sense this work may be considered as a further development of Bulaevskii et al. 
(1977) approach . The Andreev bound states appearing near the spin-active interface within the superconducting 
gap are tunable with the magnetic properties of the interface. This can result to the switch of the junction from 
"0"— to "tt"— state with changing the transmission characteristics of the interface. This approach was also applied by 
Andersson, Cuevas and Fogelstrom (2002) to study the coupling of two superconductors through a ferromagnetic dot. 
They demonstrated that the realization of the "7r" —junction is possible in this case as well. In the framework of the 
Bogoliubov-de Gennes approach Tanaka and Kashiwaya (1997) analyzed the system consisting of two superconductors 
separated by 5— functional barrier with the spin-orientation dependent height. 

Similarly to the case of S/F multilayers we may discuss the question of the existence of the S/F/S junction with 
arbitrary equilibrium phase difference tp . Naturally, the form Eq. (49) for the energy of the junction may give the 
minima at ip = and ip = tt only. A more general expression for the Josephson junction energy takes into account 
the higher order terms over the critical current which leads to the appearance of the higher harmonics over ip in the 
current-phase relationship. Up to the second harmonic, the energy is 

s = ^ (1 - cos ^2^T cos2 ^ (61) 
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and the current is 



j (ip) — I c sin ip + h sin 2ip. (62) 

If the sign of the second harmonic term is negative I2 < 0, then the transition from "0"— to "71-"— phase will be 
continuous, and the realization of the "tpo" —junction becomes possible. In general, the "</?o" — junction may exist 
if j (ifo) = and (dj/dip) Vg > 0. The calculations of the current-phase relationships for different types of S/F/S 
junctions (Golubov et al, 2004, Radovic et al., 2003 and Cayssol and Montambaux, 2004) show that (dj/dtp) < 0, 
and therefore the transition between "0"— and "71-"— states occurs to be discontinuous. 

The presence of the higher harmonics in the j (<p) relationship prevents the vanishing of the critical current at the 
transition from "0"— to "71"" — state. This is always the case when the transition occurs at low temperature. Theoretical 
studies of the properties of clean S/F/S junctions at T < T c (Buzdin et al., 1982, Chtchelkatchev et al., 2001, and 
Radovic et al., 2003) confirm this conclusion. 

Zyuzin and Spivak (2000) argued that the mesoscopic fluctuations of the critical current may produce the "7r/2" — 
superconducting Josephson junction. Such situation is possible when the thickness of F layer is close to 2d^ . The 
spatial variations of the thickness of F layer lead to the appearance of the second harmonic term in Eq. (62) with 
I2 < (Buzdin and Koshelev, 2003), and thus the realization of the 'Vo" — junction becomes possible at 2df « 2dj . 

C. Experiments with "71-" -junctions 

The temperature dependence of the critical thickness d°j is at the origin of the observed by Ryazanov et al. (2001a) 
very specific temperature dependence of the critical current I C {T) (see Fig. 14). With decreasing temperature for 
specific thicknesses of the F layer (around 27 nm), a maximum of I c is followed by a strong decrease down to zero, 
after which I c rises again. 

This was the first unambiguous experimental confirmation of the " — 7r" transition via the critical current mea- 
surements. Ryazanov et al. (2001a) explained their results by a model with a small exchange field h <~ T c . The 
Cu^Nii-a, alloy used in their experiments has the Curie temperature ~ 20 — 30JT and this implies that the ex- 
change field must be higher 100K. In consequence, it seems more probable that the thickness of the F layer was in 
the range d c j (0) < df < df (T c ), which provides the strong non- monotonous temperature dependence of I c . Also, the 
experimental estimate of £/ ~ 10 nm is too large for expected value of the exchange field. 

Recent systematic studies of the thickness dependence of the critical current in junctions with Cu^Nii-j; alloy as a 
F layer (Ryazanov et al., 2004), have revealed very strong variation of I c with the F layer thickness. Indeed, the five 
orders change of the critical current was observed in the thickness interval (12 — 26) nm. The natural explanation 
of such a strong thickness dependence is the magnetic scattering effect which is inherent to the ferromagnetic alloys. 
The presence of rather strong magnetic scattering in Cua;Nii_ x alloy S/F/S junctions was noted also by Sellicr et 
al. (2003). The magnetic scattering strengthens the decrease of the critical current with the increase of the F layer 
thickness, and at the same time it increases the period of I c (2df) oscillations. The general expression for the I c (2df) 
dependence, taking into account the magnetic scattering is given in Appendix B, Eq. (101). The attempts to describe 
the experimental data of Ryazanov et al. (2004) on the I c (2df) dependence with the help of this expression provided 
hints on the existence of the another minimum I c (2df) at smaller F layer thickness - around 10 nm. The very recent 
experiments with the junctions with the F layer thicknesses up to 7 nm have confirmed this prediction (Ryazanov et 
at, 2005) - see Fig. 15. The existence of the first "0 — 7r" transition at 2df w 11 nm means that previously reported 
transitions in Cu x Nii_ x junctions were actually the transitions from "71-" to "0"— phase (and not as was assumed, 
from "0" to "7r"— phase). It means also that now it is possible to fabricate the "it" — junctions with a 10 4 times 
higher critical current. Note, that the first measurements (Frolov et al, 2004) of the current-phase relation in S/F/S 
junction with Cuo.47Nio.53 F layer provided no evidence of the second harmonic in j (ip) relationship at the "0" - 
"tt" transition. These measurements were performed using the junction with F layer thickness around 22 nm, i. e. 
near the second minimum on the I c (2df) dependence. The much higher critical current near the first minimum (at 
2df « 11 nm) may occur to be very helpful for a search of the second harmonic. 

The results of Ryazanov et al. (2001a) on the temperature induced crossover between 0— and tt— states were recently 
confirmed in the experiments of Sellier et al. (2003). Kontos et al. (2002) observed the damped oscillations of the 
critical current as a function of F layer thickness in Nb/ Al / AI2O3/ PdNi/N b junctions. The measured critical current 
with the theoretical fit (Buzdin and Baladie, 2003) are presented in Fig. 16. Blum et al. (2002) reported the strong 
oscillations of the critical current with the F layer thickness in Nb/Cu/Ni/Cu/Nb junctions. 
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Bulaevskii et al. (1977) pointed out that "71"" -junction incorporated into a superconducting ring would generate 
a spontaneous current and a corresponding magnetic flux would be half a flux quantum $o- The appearance of the 
spontaneous current is related to the fact that the ground state of the "7r" -junction corresponds to the phase difference 
7r and so, this phase difference will generate a supercurrent in the ring which short circuits the junction. Naturally 
the spontaneous current is generated if there are any odd number of "7r" -junctions in the ring. This circumstance has 
been exploited in a elegant way by Ryazanov et al. (2001c) to provide unambiguous proof of the "7r"-phase transition. 
The authors (Ryazanov et al., 2001c) observed the half-period shift of the external magnetic field dependence of the 
transport critical current in triangular S/F/S arrays. The thickness of F layers of the S/F/S junctions was chosen 
in such a way that at high temperature the junctions were the usual "0" -junctions, and they transformed into the 
"7r" -junctions with the decrease of the temperature (Ryazanov et al., 2001a). 

Guichard et al. (2003) performed similar phase sensitive experiments using dc SQUID with "7r"-junction. The 
total current I flowing trough the SQUID is the sum of the currents I a and h flowing through the two junctions, 
I = I a +Ib- If the junctions have the same critical currents I c and both are "0" -junctions, then I a = I c sunp a and 
h = Ic sin (fib, where tp a and tpi, are the phase differences across the junctions. Neglecting the inductance of the loop 
of SQUID, the phase differences satisfy the usual relation (Barone and Paterno, 1982), ip a — <Pb = 27r$/$ 7 w here $ is 
the flux of the external magnetic field through the loop of the SQUID. The maximum critical current of the SQUID 
will be I max = 21 c cos(7T$/$ )- In the case when one of the junctions (let us say b) is the "7r" -junction with the 
same critical current, the current flowing through it /(, = — I c sin ipb = I c sin (ipb + i") ■ Therefore the maximum critical 
current of the SQUID in this case will be i^ax = 2/ c cos(7r$/$o + 1"/2), and the diffraction pattern will be shifted of 
half a quantum flux. If both junctions are the "7r" -junctions the diffraction pattern will be identical to the diffraction 
pattern of the SQUID with two "0" -junctions. Namely this was observed on experiment by Guichard et al. (2003) 
with SQUID containing junctions with PdNi ferromagnetic layers, see Fig. 17. 

Recently Bauer et al. (2004) measured with the help of micro Hall-sensor the magnetization of a mesoscopic 
superconducting loop containing a PdNi ferromagnetic 'V -junction. These measurements also provided a direct 
evidence of the spontaneous current induced by the 'V -junction. 

VI. COMPLEX S/F STRUCTURES 

A. F/S/F spin-valve sandwiches 

The strong proximity effect in superconductor-metallic ferromagnet structures could lead to the phenomenon of 
spin-orientation-dependent superconductivity in F/S/F spin-valve sandwiches. Such type of behavior was predicted 
by Buzdin et al. (1999) and Tagirov (1999) and recently has been observed on experiment by Gu et al. (2002). 
Note that a long time ago De Gennes (1966b) considered theoretically the system consisting of a thin S layer in 
between two ferromagnetic insulators. He argued that the parallel orientation of the magnetic moments is more 
harmful for superconductivity because of the presence of the non-zero averaged exchange field acting on the surface 
of the superconductor. This prediction has been confirmed on experiment by Hauscr (1969) on In film sandwiched 
between two Fe304 films and Deutscher and Meunier (1969), on a In film between oxidized FeNi and Ni layers, see 
Fig. 18. Curiously, the experiments of Deutscher and Meunier (1969) correspond more to the case of the metallic 
F/S/F sandwiches as the authors report rather low interface resistance. 

To consider the spin-orientation effect in metallic F/S/F sandwiches we use the notations analogous to that of 
section 4. More precisely, to have a direct connection with the corresponding formula of Section 4, we assume that 
the thickness of the F layers is df and the S layer - 2d s , see Fig. 19. 

Also, to provide a simple theoretical description we consider the case d s -C £ s with only two orientations of the 
ferromagnetic moments: parallel and antiparallel. The case of arbitrary orientations of the ferromagnetic moments 
needs the introduction of triplet components of the anomalous Green's functions. The first attempt of such analysis 
was made by Baladie et al. (2001), but on the basis of the incomplete form of the Usadel equation. The full correct 
calculations for this case has been performed by Volkov et al. (2003), Bergeret et al. (2003), and Fominov et al. 
(2003a). 

In fact, we only need to analyze the case of the antiparallel orientation of the ferromagnetic moments because the 
case of the parallel orientation is completely equivalent to the "0"— phase in S/F multilayered structure (Section 4) 
with the F layers two times thinner than in a F/S/F sandwich. In other words, our choice of notations permits for the 
parallel orientation case to use directly the corresponding expressions for the critical temperature for the "0"— phase 
from Section 4. To analyze the antiparallel orientation case, we follow the approach used in Section 4, but we need 
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to keep the linear over x term in the expansion of the anomalous Green's function in the S layer Eq. (35) 

ft, 2 



F s (x,u>) = F ( 1 + a u x - —x 



(63) 



With the help of the Usadel equation (31), we readily find that Fo has the form (36) with the pair-breaking parameter 
t.7 1 determined by the expression 
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Let us suppose that the exchange field is positive (+h) in the right F layer and then for d s + df > x > d s 



Ff(x,u > 0) = ,4 cosh 







(x-d s - d f ) 



while for the left F layer, the exchange field is negative and for — d s — df < x < — d s we have 



Ff(x,u> > 0) = Bcosh 



(x + d s + df) 



(64) 
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Taking into account the explicit form of the function Ff(x) and the boundary conditions (32), we may see that for the 

and the pair-breaking parameter for this case r s _1 = t~\ p may be 
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The second term in the right-hand side of the eq. (67) may be important only in the limit of very small df and 
we will omit it further. The boundary conditions Eqs. (32) permit us to calculate the parameter r s _1 , provided the 
anomalous Green function in the F layer is known. For the parallel alignment of the ferromagnetic moments it is just 
r~p = t\To , where t~q is given by the Eq. (40), while for the antiparallel alignment it is just 
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(68) 



In result, we obtain the following simple formula for the critical temperature T c p for the parallel orientation and T c 
for the antiparallel one 
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(69) 



in ^=K0-*{^ +Re (^^)}- (70) 

The different kinds of T c (df) curves are presented in Fig. 20. 

We see that the interface transparency is the important factor, controlling the spin- valve effect in F/S/F structures. 
It is interesting that the optimum condition for the observation of this effect in the case of the non-negligeable interface 
transparency is the choice df ~ (0.1 — 0.4) £/. 

In the case when the F layer thickness exceeds £/, the critical temperature practically does not depend on df. This 
case for the transparent S/F interface (-f B = 0) was considered by Buzdin et al. (1999), and the critical temperatures 
for the parallel and antiparallel aligncmcnts are presented in Fig. 21. The finite interface transparency strongly 
decreases the spin- valve effect, and for the parameter 7 y B > 5 the dependence of the critical temperatures on the 
mutual orientation of ferromagnetic moments is hardly observable. 
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The thermodynamic characteristics of F/S/F systems were studied theoretically by Baladie and Buzdin (2003) 
and Tollis (2004) in the framework of Usadel formalism and it was noted that the superconductivity always remains 
gapless. 

Bagrets et al. (2003) developed a microscopic theory of F/S/F systems based on the direct solution of the Gor'kov 
equations for the normal and anomalous Green's functions. The main mechanism of the electron scattering in F 
layers was supposed to be of the s — d type. The results of this microscopical analysis were in accordance with the 
quasiclassical approach and provided a reasonable quantitative description of the experimental data of Obi et al. 
(1999) on T c (df) dependence in Nb/Co multilayers. 

Krunavakarn et al. (2004) generalized the approach of Fominov et al. (2002) to perform exact numerical calculations 
of the nonmonotonic critical temperature in F/S/F sandwiches. They demonstrated also that the Takahashi-Tachiki 
(1986) theory of the proximity effect is equivalent to the approach based on the Usadel equations. 

Bozovic and Radovic (2002) studied theoretically the coherent transport current through F/S/F double-barrier 
junctions. The exchange field and the interface barrier reduce the Andreev reflection due to the enhancement of the 
normal reflection. Interestingly, that the conductance is always higher for parallel alignment of the ferromagnetic 
moments. The similar conclusion was obtained in work of Yamashita et al. (2003). Such behavior is related with the 
larger transmission for the normal tunneling current in this orientation. The calculations also revealed the periodic 
vanishing of Andreev reflection at the energies of geometrical resonance above the superconducting gap. 

The case of insulating F layers (De Gennes, 1966b) corresponds to the situation when the superconducting electrons 
feel the exchange field only on the surface of S layer. We may describe this case taking formally the limit df — > 
with TgQ 1 = ih-j-, where a is the distance of the order of the interatomic one, which describes the region near 
the S/F interface where the exchange interaction (described by the exchange field h ) with electron spins takes 
place. In fact it simply means that, for the parallel orientation case, the superconductor is under the influence of 
the averaged exchange field h = h-^-, while for the antiparallel orientation this field is absent. Careful theoretical 
analysis of the system consisting of the superconducting film sandwiched between two ferromagnetic semiconducting 
insulators with differently oriented magnetization was performed by Kulic and Endres (2000) for both singlet and 
triplet superconductivity cases. In the case of a triplet superconductivity, the critical temperature depends not only 
on the relative orientation of the magnetization but also on its absolute orientation. 



B. S-F-l-F'-S heterostructures and triplet proximity effect 

A bunch of theoretical works was devoted to the analysis of more complex S/F systems. Proshin et al. (2001) 
(see also Izyumov et al. 2002) studied the critical temperature of S/F multilayers with alternating magnetization of 
adjacent F layers. The same authors (Izyumov et al, 2000 and Izyumov et al. (2002)) also proposed the 3D LOFF 
state in F/S contacts. However, this conclusion was based on controversial boundary conditions, corresponding to the 
different in plane 2D wave-vectors on the both sides of the contact - see the comment by Fominov et al. (2003b) and 
the reply of Khusainov and Proshin (2003). 

Koshina and Krivoruchko (2001) (see also Golubov et al. 2002a) studied the Josephson current of two proximity 
S/F bilayers separated by an insulating (I) barrier and demonstrated that in such S/F-I-F/S contact the 7r-phase may 
appear even at very small F layer thickness (smaller than £/). The mechanism of the 7r-phase transition in this case 
is related to the rotation on n/2 of the phase of the anomalous Green's function F on the S/F boundary in addition 
to the jump of its modulus. To demonstrate this we consider the thin F layer of the thickness df << £ s in contact 
with a superconductor. If the x = corresponds to the S/F interface, and x — df is the outer surface of the F layer, 
then the solution of the linearized Usadel equation in the ferromagnet is 



Ff (x,u> > 0) = A cosh 
Using the boundary condition Eq.(32) we may easily obtain 
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(72) 



In the case of a rather low interface transparency, 7sCn^//C/ >> 1> the jump of the phase of the F function at the 
interface is practically equal to — it/2 : 
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Ff (0, u > 0) w F s (0, w > 0) exp(-i-) 




(73) 



Koshina and Krivoruchko (2001) and Golubov et al. (2002a) argued that at each S/F interface in the S/F-I-F/S 
contact the phase jump — ir/2 occurs, and the total phase jump in the equilibrium state would be ir. 

Kulic and Kulic (2001) calculated the Joscphson current between two superconductors with a helicoidal magnetic 
structure. They found that the critical current depends on the simple manner on the relative orientation 9 of the 
magnetic moments on the banks of contact : 



where R_(R + ) corresponds to the same (opposite) helicity of the magnetization in the banks. Depending on the 
parameters of the helicoidal ordering, the value of R± may be either smaller or larger than 1. If R± > 1, than I c may 
be negative for some misorientation angles 9, which means the realization of the 7r - phase. Interestingly that tuning 
the magnetic phase 9, it is possible to provoke a switch between - and tt - phase. As it may be seen from Eq. (74), 
the critical current of the Josephson junction is maximal for the antiparallel orientation (9 — tt) of the magnetizations 
in the banks. 

Bergeret et al. (2001a) studied the Joscphson current between two S/F bilayers and pointed out the enhancement 
of the critical current for an antiparallel alignment of the ferromagnetic moments. They demonstrated that at low 
temperatures the critical current in a S/F-I-F/S junction may become even larger than in the absence of the exchange 
field (i. e. if the ferromagnetic layers are replaced by the normal metal layers with h = 0). More in details (taking into 
account different transparency of S/F interfaces and different orientations of the magnetization in the banks) these 
junctions were studied theoretically by Krivoruchko and Koshina (2001), Golubov et al. (2002a), Chtchelkatchcv et 
al. (2002) and Li et al. (2002). Blanter and Hekking (2004) used Eilenberger and Usadel equations to calculate the 
current-phase relation of Josephson junction with the composite F layer, consisting of two ferromagnets with opposite 
magnetizations. 

Bergeret et al. (2001b) and Kadigrobov et al. (2001) analyzed in the framework of Usadel equations the proximity 
effect in S/F structures with local inhomogeneity of the magnetization. They obtained an interesting conclusion that 
the varying in space magnetization generates the triplet component of the anomalous Green's function (~ (\I r f\f r f)) 
which may penetrate in the ferromagnet at distances much larger than It is not however the triplet superconduc- 
tivity itself because the corresponding triplet order parameter would be equal to zero, unlike the superfluidity in He 3 , 
for example. In general, the triplet components of the anomalous Green's function always appear at the description of 
the singlet superconductivity in the presence of rotating in space exchange field. For example, they were introduced 
by Bulaevskii et al. (1980) in the theory of coexistence of superconductivity with helicoidal magnetic order. An 
important finding of Bergeret et al. (2001b) and Kadigrobov et al. (2001) was the demonstration of the fact that in 
some sense the triplet component is insensitive to the pair-breaking by the exchange field. Therefore its characteristic 



decaying length is the same as in the normal metal, i. e. = V^r- The triplet long-range proximity effect 



could explain the experiments on S/F mesoscopic structures (Giroud et al, 1998 and Petrashov et al., 1999), where 
a considerable increase of the conductunce below the superconducting critical temperature was observed at distances 
much larger than £/. 

In their subsequent works Bergeret et al. (2003) and Volkov et al. (2003) studied the unusual manifestation of this 
triplet component in S/F multilayered structures. The most striking effect is the peculiar dependence of the critical 
current in multilayered S/F structures on the relative orientation of the ferromagnetic moments. For the collinear 
orientation, the triplet component is absent, and provided the thickness of the ferromagnetic layer df >> , the 
critical current is exponentially small. On the other hand, if the orientation of the magnetic moments is noncollincar 
then the triplet component of the superconducting condensate appears. Its decaying length ^x,d is much larger than 
£/, and namely this triplet component realizes the coupling between the adjacent superconducting layers. When the 
thicknesses of F layers are in the interval of £-r,d >> df » £/, then this coupling occurs to be strong. In result, 
the critical current is maximal for the perpendicular orientation of the adjacent ferromagnetic moments, and it may 
exceed many times the critical current for their parallel orientation. Due to the mesoscopic fluctuations (Zyuzin 
et al. 2003), the decay of the critical current for collinear orientation of the magnetic moments is not exponential. 
Nevertheless, for this orientation it would be very small, and this circumstance do not change the main conclusion on 
the existence of the long range triplet proximity effect. A lot of interesting physics is expected to emerge in the case 



Ic = IcO (1 - R± cosfl) , 



(74) 




25 



of S/F systems with genuine triplet superconductors. For example, the proximity effect would be strongly dependent 
on the mutual orientation of the magnetic moments of the Cooper pairs and ferromagnets. 

The long range triplet proximity effect was predicted to exist in the dirty limit. An interesting question is how it 
evolves in the clean limit. In this regime there is no characteristic decaying length for the anomalous Green's function 
in a ferromagnet (see Eqs. (21), (22)), and the angular behavior of the critical current in S/F multilayers may be quite 
different. If, for example, we apply the Eilcnberger equations for the description of clean S/F/F'/S structure with 
antiparallel ferromagnetic layers with equal thicknesses, the exchange field completely drops (Blanter and Hckking, 
2004). Therefore, the critical current will be the same as for the non magnetic interlayers. In this case it is difficult 
to believe that for the perpendicular orientation of the magnetic moments the critical current could be even higher. 
The microscopical calculations in the framework of the Bogoliubov-de Gennes equations of the properties of S/F 
multilayers with non-collincar orientation of the magnetic moments would be of substantial interest. 

Barash et al (2002) studied the Josephson current in S-FIF-S junctions in clean limit within the quasiclassical theory 
of superconductivity, based on the so-called Ricatti parametrization (Schopol and Maki, 1995). They obtained the 
striking nonmonotonic dependences of the critical current on the misorientation angle of the ferromagnetic moments. 
However, even for a rather high transparency of I barrier (D = 0.8), the maximum of the critical current occurred for 
the antiparallel orientation of the magnetic moments. 

VII. ATOMIC THICKNESS S/F MULTILAYERS 

A. Layered ferromagnetic superconductors 

In this section, we consider an atomic-scale multilayer F/S system, where the superconducting (S) and the ferromag- 
netic (F) layers alternate. When the electron transfer integral between the S and F layers is small, superconductivity 
can coexist with ferromagnetism in the adjacent layers. Andreev et al, (1991) demonstrated that the exchange field 
in F layers favors the 'V — phase behavior of superconductivity, when the superconducting order parameter alternates 
its sign on the adjacent S layers. 

Nowdays several type of layered compounds, where superconducting and magnetic layers alternate, are known. 
For example in Smi.85Ceo.i5Cu04 (Sumarlin et al, 1992), which reveals superconductivity at T c = 23.5 K, the 
superconducting layers are separated by two ferromagnetic layers with opposite orientations of the magnetic moments 
and the Neel temperature is T/v = 5.9 K. Several years ago, a new class of magnetic superconductors based on the 
layered perovskite ruthenocuprate compound RuSr 2 GdCu20s comprising Cu02 bilayers and Ru02 monolayers has 
been syntesized (see for example McLaughlin et al., 1999 and references cited there). In RuS^GdC^Os, the magnetic 
transition occurs at Tm ~ 130 — 140, K and superconductivity appears at T c ~ 30 — 50 K. Recent measurements of the 
interlayer current in small-sized RuSr 2 GdCus single crystals showed the intrinsic Josephson effect (Nachtrab et al, 
2004). Apparently, it is a week ferromagnetic order which is realized in this compound. Though the magnetization 
measurements give evidence of the small ferromagnetic component, the neutron diffraction data on RuS^GdCu^Os 
(Lynn et al, 1992) revealed the dominant antifcrromagnetic ordering in all three directions. Later, the presence 
of ferromagnetic in-plane component of about (0.1-0.3)/xs has been confirmed by neutron scattering on isostructural 
RuS^YCu^Os (Tokunaga et al., 2001). In addition, in the external magnetic field the ferromagnetic component grows 
rapidly at the expense of the antiferromagnetic one. 

Due to the progress of methods of the multilayer preparation, the fabrication of artificial atomic-scale S/F superlat- 
tices becomes possible. An important example is the YBa2Cu307/La 2 /3Ca 1 / 3 Mn03 superlattices (Sefrioui et al, 2003 
and Holden et al, 2003). The manganite half metallic compound La2/3Cai/ 3 Mn03 (LCMO) exhibits colossal magne- 
toresistance and its Curie temperature 9 = 240 K . The cuprate high-T c superconductor YBa2Cu307 (YBaCuO) with 
T c = 92 K, have the similar lattice constant as LCMO which permits to prepare the very high quality YBaCuO/LCMO 
superlattices with different ratio of F and S layers thicknesses. The proximity effect in these superlattices occurs to 
be extremely long-ranged. For a fixed thickness of the superconducting layer, the critical temperature is dependent 
over a thickness of LCMO layer in the 100 nm range (Sefrioui et al, 2003 and Pena et al, 2004). This is very unusual 
behavior because the YBaCuO and LCMO are strongly anisotropic layered systems with very small coherence length 
in the direction perpendicular to the layers (0.1 — 0.3 nm). Somewhat similar giant proximity effect has been recently 
reported in the non-magnetic trilayer junctions Lai.85Sro.i5Cu04/La2Cu04+d/Lai.85Sro.i5Cu04 (Bozovic et al, 2004) 
and in the superconductor-antiferromagnet YBa 2 Cu 3 07/ Lao.4sCao.55Mn03 superlattices (Pang et al, 2004). The 
observed giant proximity effect defies the conventional explanations. Bozovic et al. (2004) suggested that it may be 
related with resonant tunneling, but at the moment the question about the nature of this effect is open. 
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B. Exactly solvable model of the 'V'-phase 



Let us consider the exactly solvable model (Andreev et al, 1991) of alternating superconducting and ferromagnetic 
atomic metallic layers. For simplicity, we assume that the electron's motion inside the F and S layers is described by 
the same energy spectrum £ (p). Three basic parameters characterize the system : t is the transfer energy between 
the F and S layers, A is the Cooper pairing constant which is assumed to be non zero in S layers only, and h is the 
constant exchange field in the F layers. The Hamiltonian of the system can be written as 

H= ^2 £(p)a+ i<T (p)<w(p) + Hi n ti +H int2 + 

+* [«m ff (PH,-^(P) + <+l,- v (P)am ff (p) + h.c] , 

H intl = (75) 

| a nl CT (Pl) a n,l,-<T(-Pl) a ™a,- CT (-P2)ar l l ( T(P2), 
Pi,P2,n,cr 

H lnt2 =~h^2 aa n ,-l,a(P) a n ,-1,<t(p), 
p,n,cr 

where a^ i(7 is the creation operator of an electron with spin a in the n th elementary cell and a momentum p in the 
layer i, where i = 1 for the S layer, and i = — 1 for the F layer, and g is the pairing constant. The important advantage 
of this model is that the quasiparticle Green's functions can be calculated exactly and the complete analysis of the 
superconducting characteristic is possible. Assuming that the order parameter changes from cell to cell in the manner 
A„ = |A| e lkn , the self-consistency equation for the order parameter |A| reads 
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where A = gN(0) and ui± — iuj±^ (p) ,w± = u± + h. The quasimomentum q lies in the direction perpendicular to the 
layers, and T q = 2t cos (q/2) e lq / 2 . In the limit of a small transfer integral t « T c , where T c is the bare mean-field 
critical temperature of the S layer in the absence of coupling (t — 0), we arrive at the following equation for the 
critical temperature T* : 

In Tl = -nT^t 2 V , - t — + (77) 
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The critical temperature T* is close to the bare critical temperature T c and as is seen from Eq. (77), for h = 0, the 
maximal T* corresponds to k = 0, i.e. the superconducting order parameter is the same at all layers. It is worth 
to note that as the exchange field on the F layers grows, tunneling becomes energetically more costly, so the leading 
term second order in t falls as 1/h 2 for large h and the critical temperature increases. This is related to the fact that, 
due to the decrease of the coupling the effective exchange field induced on the S layers decreases with the increase 
of h. For h >> T c , the coefficient of the cosfc term has a negative sign and the maximal T* corresponds to k = it, 
so the transition occurs to the 7r-phase with an alternating order parameter A„ = |A| (—1)™. Numerical calculations 
(Andreev et at, 1991) give for the critical value of the exchange field (at which k changes from to n) h c = 3.77T C , 
and the complete (h,T) phase diagram is presented in Fig. 22. 

At T = the transition to the "7r"-phase occurs at h cn = 0.87T C . The analysis of Prokic et al. (1999) and Houzet 
et al. (2001) shows that the perpendicular critical current vanishes at the line of the transition from the "0"- to the 
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"7r" -phase and the Josephson coupled superconducting planes are decoupled. Strictly speaking, the critical current 
vanishes only in ~ t 4 approximation, see Eq. (77). The term ~ t 8 gives the contribution ~ t 8 cos2/c, and the critical 
current at the transition to the 'V -phase will drop to the very small value ~ I c (t/T c ) 8 . Note that the sign of the 
second harmonic in j((p) relation generated by this ~ t 8 term is positive, and therefore the transition from "0"- to 
the "7r" -phase is discontinuous. 

In result, if the exchange field is in the interval h c o < h < 3.77T C , the "0-7r" transition may be easily observed 
with the lowering of the temperature due to the nonmonotoneous behavior of the Josephson plasma frequency and 
the parallel London penetration (Houzet et al, 2001). However the typical value of the exchange field is rather high 
and more probable is the situation h » T Cl and so the system will be in the 'V-phasc at any temperatures. This 
is consistent with the recent experiments of Nachtrab et al. (2004) on RuS^GdC^Os presenting no evidence of 
superconducting planes decoupling with temperature. In RuS^GdC^Os, the superconducting pairing is probably of 
the d-wave type. This case was analyzed theoretically by Prokic and Dobrosavljevic-Grujic (1999), and the scenario of 
the "7r" -phase appearence is very close to the case of the s-wave superconductivity. Calculations of electronic density 
of states by Prokic and Dobrosavljevic-Grujic (1999) and Prokic et al. (1999) revealed some changes inherent to the " 
0-7r" transition, but, apparently, the experimental identification of the 7r-phase in the atomic-scale S/F superlattices 
is an extremely difficult task. In principle, if the superlattice consists of an even number of superconducting layers, 
then the phase of the order parameter at the ends will differ by 7r, and the entire system will function as a Josephson 
"7r" -junction. The spontaneous current in a superconducting loop containing such a 'V-junction could be observed 
at an experiment analogous to the one made by Bauer et al. (2004). 

The model Eq.(75) permits to analyze the transition from the quasi-2D to 3D system with the increase of the 
transfer intergral t. At t < T c , instead of the "7r" -phase, the LOFF state with modulation along the superconducting 
layers appears and the system becomes analogous to the 3D superconductor in an uniform exchange field (Houzet 
and Buzdin, 2002). 

Buzdin and Daumens (2003) considered the spin waive effect in the F/S/F structure consisting of three atomic layers 
and described by the model Eq. (75). Analogously to the F/S/F spin- waive sandwiches (see Section 6), the critical 
temperature is maximal for the antiparallel orientation of the ferromagnetic moments. However, at low temperature, 
the situation is inversed. Namely, the superconducting gap occurs to be larger for the parallel orientation of the 
ferromagnetic moments. This counter-intuitive result of the inversion of the proximity effect may be understood on 
the example of the ferromagnetic half-metal. Indeed at T = 0, the disappearance of the Cooper pair in a S layer 
means that two electrons with opposite spin must leave it. If the neighbouring F layers of half- metals are parallel, 
then, for one spin orientation, they are both insulators and the electron with this spin orientation can not enter it. 
It results in the impossibility of the pair destruction. On the other hand, for the antiparallel orientation of the F 
layers, for any electron spin orientation there is an ajacent normal layer and a Cooper pair can leave the S layer. 
Such behavior contrasts with the diffusive model prediction (Baladie and Buzdin, 2003 and Tollis, 2004) but is in 
accordance with the T — results obtained in the framework of the multiterminal model for S/F hybrid structures 
(Apinyan and Melin, 2002). Apparently, it is a special property of the clean limit of the atomic-layer S/F model, and 
it disappears in the case of several consequitive S layers per unit cell (Melin and Feinberg, 2004) . 

VIII. SUPERCONDUCTIVITY NEAR THE DOMAIN WALL 

In the previous discussion of the properties of S/F heterostructures, we have implicitly assumed that the ferromagnet 
has uniform magnetization, i. e. there are no domains. It practice, the domains appear in ferromagnets quite easily 
and special conditions are usually needed to obtain the monodomain ferromagnet. In standard situation, the size 
of the domains is much larger than the superconducting coherence length, and << £ s , therefore the Cooper pair 
will sample the uniform exchange field. However, a special situation with the S/F proximity effect is realized near 
the domain wall, where the magnetic moments and the exchange field rotate. The Cooper pairs feel the exchange 
field averaged over the superconducting coherence length. Naturally, such averaged field will be smaller near the 
domain wall, which leads to the local decrease of the pair-breaking parameter. As the result, we may expect that 
superconductivity would be more robust near the domain wall. In particular, the critical temperature T cw for the 
superconductivity localized near the domain wall would be higher than that of the uniform S/F bilayer T*. For 
bulk ferromagnetic superconductors, the critical temperature of the superconductivity localized near the domain wall 
was calculated by Buzdin et al, (1984). The experimental manifestations of the domain wall superconductivity in 
Nio.80Feo.20/Nb bilayers (with Nb thickness around 20 nm) were observed by Rusanov et al. (2004). The Neel-type 
domain walls in Permalloy (Ni .8oFeo.2o) are responsible for the local increase of the critical temperature around 10 
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mK. The width of the domain walls w in Permalloy films used in (Rusanov et al., 2004) is rather large w ~ 0.5 fim, 
i. e. much larger than the superconducting coherence length of niobium. The rotation angle a of the exchange field 
at the distance £ s may be estimated asa~( s /t() , and so the averaged exchange field h av is slightly smaller than the 
field h far away from the domain wall: (h — h av ) /h ~ (^ s /w) 2 . Therefore, the relative decrease of the pair-breaking 
parameter 1-7 1 in Eq. (40) will be also of the order ~ (^ s /w) . From Eqs. (40, 43) we obtain the following estimate 
of the local increase of the critical temperature 

^V^~(6M 2 , (78) 

1 c 

which is of the same order of magnitude as the effect observed on the Nio.80Feo.20/Nb bilayers. Keeping in mind the 
temperature dependence of the superconducting coherence length £(T) ~ ^y ^j^ -Tj, we see that the condition of 

the domain wall superconductivity appearance is simply £,(T CW ) ~ w. 

In the case of a very thin domain wall, the variation of the exchange field is a step-like and the local suppression of 
the pair-breaking parameter occurs at the small distance of the order £f << £ s near the domain wall. The situation 
resembles the enhancement of the superconducting pairing near the twin planes (Khlyustikov and Buzdin, 1987). The 
variation of the pair-breaking occuring over a distance £/ induces a superconducting order parameter over a distance 
£.(T CW ) near the domain wall and the effective relative decrease of the pair-breaking parameter will be of the order 
of £f/£(T cw ). Therefore, if the shift of the critical temperature of the S/F bilayer is comparable with T c itself, i. 
c. (T c — T*) /T c ~ 1, the critical temperature T cw of the superconductivity, localized near the domain wall may be 
estimated from the condition Tc ^ T , Tc ~ £f/£(T cw ). In result we have 

~ (0V6) 2 , (79) 



which is around (1-5)% for typical values of £/ and £ s . A small width of the domain walls is expected in experiments 
of Kinscy, Burnell, and Blamirc (2001) on the critical current measurements of Nb/Co bilayers. The domain walls 
occured to be responsible for the critical current enhancement below T* = (5. 24 ±0.05) K. In the presence of domains 
walls the non-zero critical current has been observed at (5.4 ± 0.05) K, slightly above T* . 

It is worth to note that the effect of the increase of the critical temperature in the vicinity of a domain wall is weak 
for very large and very thin domain wall. The optimum thickness, when the effect may be ralatively strong is w ~ £ s . 

In the case of a perpendicular easy-axis the branching of the domains may occur near the surface of magnetic 
film. If the scale of this branching is smaller than the superconducting coherence length, the effective exchange field 
is averaged, and the pair breaking parameter will be strongly decreased. This mechanism has been proposed in 
(Buzdin, 1985) to explain the presence of traces of superconductivity at low temperature in re-entrant ferromagnetic 
superconductors. The similar effect may take place in S/F bilayers and in such case the superconductivity would be 
extremely sensitive to the domain structure. Rather weak magnetic field would suffice to modify the branching of 
domains and supress superconductivity. 

Up to now we have concentrated on the interplay between superconductivity and ferromagnetism caused by the 
proximity effect related to the passing of electrons across the S/F interface. However, if the magnetic field created 
by the ferromagnet penetrates into a superconductor, it switches on the orbital mechanism of superconductivity 
and magnetism interaction. The situation when it is the only one mechanism of superconductivity and magnetism 
interaction is naturally realized in the case, when the ferromagnet is an insulator, or the buffer oxide layer separates 
the superconductor and the ferromagnet. The hybrid S/F systems have been intensively studied in connection with 
the problem of the controlled flux pinning. Enhancement of the critical current has been observed experimentally 
for superconducting films with arrays of submicron magnetic dots and antidots(see, for example Van Bael et al., 
2002a and Van Bael et al, 2002b, and references cited therein), and for S/F bilayers with a domain structure in 
ferromagnetic films (Garcia-Santiago et al, 2000). A theory of vortex structures and pinning in S/F systems at rather 
low magnetic field has been elaborated by Lyuksyutov and Pokrovsky (1998), Bulaevskii et al. (2000), Erdin et al. 
(2002) and Milosevic et al., (2002a). This subject is discussed in details in the recent review by Lyuksyutov and 
Pokrovsky (2004). 

The nucleation of the superconductivity in the presence of domain structure has been theoretically studied by 
Buzdin and Melnikov (2003), and Aladyshkin et al. (2003) in the case of magnetic film with perpendicular anisotropy. 
The conditions of the superconductivity appearance occur to be more favorable near the domain walls. Recently the 
manifistation of the domain wall superconductivity was revealed on experiment by Yang et al. (2004). They deposited 
on the single crystal ferromagnetic BaFe^Oig substrate a 10 nm Si buffer layer and then a 50 nm Nb film. The 
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strong magnetic anisotropy of BaFe^Oig assures that its magnetisation is perpendicular to the Nb film. The very 
characteristic R(T) dependences and pronounced hysteresis effects have been found in the resistance measurements 
in the applied field. 

A different situation is realized if the magnetization of F layer is lying in the plane (parallel magnetic anisotropy). 
Then any type of the domain walls will be a source of the magnetic field for the adjacent S layer, and the domain wall 
locally weakens superconductivity. This idea was proposed by Sonin (1988) to create in a S layer a superconducting 
weak link (Josephson junction) attached to the domain wall. 

Lange et al. (2003) used a nanoengineered lattice of magnetic dots on the top of the superconducting film for the 
observation of the field-induced superconductivity. The applied external magnetic field provided the compensation 
of the magnetic field of the dots and increased the critical temperature. The idea of such compensation effect was 
proposed a long time ago by Ginzburg (1956) for the case of the ferromagnetic superconductors. 

The analysis of the superconducting states appearing near the magnetic dots (when the upper critical field depends 
on the angular momentum of the superconducting nucleus wave function) was done in the works of Cheng and Fertig 
(1999) and Milosevic et al. (2002b). 

IX. MODIFICATION OF FERROMAGNETIC ORDER BY SUPERCONDUCTIVITY 

A. Effective exchange field in thin S/F bilayers 

The influence of ferromagnetism on superconductivity is strong, and it leads to many experimentally observed 
consequences. Wherther the inverse is true also ? In other words, can superconductivity affect or even destroy ferro- 
magnetism ? To address this question, we start with comparing the characteristic energy scales for superconducting 
and magnetic transitions. The energy gain per atom at the magnetic transition is of the order of the Curie tem- 
perature 0. On the other hand the condensation energy per electron at the superconducting transition (Eq. (2)) is 
much smaller than T c , and it is only about ~ T c (T c /Ep) « T c . Usually the Curie temperature is higher than T c 
and ferromagnetism occurs to be much more robust compared with superconductivity. Therefore the superconduc- 
tivity can hardly destroy the ferromagnetism, but it may nevertheless modify it, if such modification do not cost too 
much energy. The example is the bulk ferromagnetic superconductors ErRh 4 B 4 , HoMoeSs and HoMogSes, where, 
in superconducting phase, ferromagnetism is transformed into a domain phase with the domain size smaller than 
the superconducting coherence length £ s (Maple and Fisher, 1982; Bulaevskii et al, 1985). Similar effect has been 
predicted by Buzdin and Bulaevskii (1988) for a thin ferromagnetic film on the surface of a superconductor. To 
illustrate this effect, we consider the S/F bilayer with S layer thickness d s smaller than the superconducting coherence 
length £ s and the F layer thickness df « << d Sl see Fig. 23. 

In the case of a transparent S/F interface, the pair-breaking parameter is given by the Eq. (40), and it is 



which simply means that the effective exchange field in the superconductor h « h-j- yjf^-^j ■ The condition of a 
transparent interface implies that the Fermi momenta are equals in both materials and this permits us to write the 
effective field as 



where vf s and VFf are the Fermi velocities in S and F layers respectively. Note however that for strong ferromagnets 
the condition of perfect transparency of the interface is different, vf^vfi — v^, where vf] and vfi are the Fermi 
velocities for two spin polarizations in ferromagnet (Zutic and Vails, 1999, and Zutic et al, 2004). 

In fact, in the considered casejaf thin F and S layers the situation is analogous to the magnetic superconductors 
with an effective exchange field h, which may also depend on the coordinates {y, z) in the plane of bilayer. Let us 
demonstrate this important point. Keeping in mind the domain structure, (see Fig. 23), where the exchange field 
depends only on the z coordinate, we may write the Usadel equations in F and S layers 



L>f d s a s 



(80) 




h = h (d f /d s ) (vFs/vFf) 
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Now let us perform the averaging procedure by integrating these equations over x. Due to the small thicknesses of F 
and S layers, the Green's functions G and F vary little with x and may be considered as constants. The integration 
of the terms with the second derivatives on x will generate Of- and Op- terms taken at the interfaces. At the interfaces 

to Ox Ox 

with vacuum these derivatives vanish and the boundary conditions Eq.(32) permit us to rely on the derivatives of F 
function on both sides of the S/F interface (the same relation is true for the G function, due to the normalization 
condition Eq. (98)). Excluding the derivatives and (^f) d , we obtain the standard Usadel equation but for 

the averaged (over the S layer thickness) Green's functions F and G 
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where the effective field h(z) = = ^T"ffj an d the condition df/d s << 1 is used to neglect the small 

renormalization of D a and u>. The possibility to introduce the effective field h(z) in the case of a thin bilaycr is quite 
natural and rather general. The same effective field may be introduced in the framework of Eilcnbcrgcr equations. 



B. Domain structure 

In the case of the uniform ferromagnetic ordering in the F layer, superconductivity can exist only if h does not 
exceed the paramagnetic limit: h < 1.24T C . This means that the thickness of the F layer must be extremely small 
df < (T c /h) d s ; even for d s ~ £ s , taking T c ~ 10 K and h ~ 5000 K, the maximum thickness of F layer only around 
1 nm. However, the ferromagnetic superconductors (Maple and Fisher, 1982; Bulaevskii et al., 1985) give us the 
example of domain coexistence phases with the exchange field larger than the paramagnetic limit. 

We may apply the theory of magnetic superconductors (Bulaevskii et al, 1985) to the description of the domain 
structure with wave vector Q >> t^^T 1 in the S/F bilayer, Fig. 23. The pair-breaking parameter associated with the 
domain structure is t^ 1 ~ ^ (Bulaevskii et al, 1985), where v — vp s is the Fermi velocity in S layer. Let us write 
the domain wall energy per unit area as a/na 2 , where a is the interatomic distance. The domain wall energy in the 
F film per unit length of the wall will be df (cr /7ra 2 ). Note that we consider the case of relatively small domain wall 
thickness w << Q _1 <<^ s and the constant <r, describing the domain wall energy is of the order of Curie temperature 
6 for the atomic thickness domain wall but may be smaller for the thick domain wall. The change of the density 
of the superconducting condensation energy due to the pair-breaking effect of domain structure is of the order of 
N(0)A 2 / (At s ). Therefore the density (per unit area) of the energy Eds related to the domain structure reads 



E DS ~N(0)d.A^ +df^T- (85) 



Its minimum is reached at 



, _ d s N(0)Aa 2 h 2 1 d s h 2 



Q 2 = n -^"" ( 86 ) 

df (tv a£o df aEp 

where £0 = Hv/ (nA). The factor which favors the existence of the domain structure is the superconducting con- 
densation energy E s ~ — N(0)d s A 2 per unit area. The domain structure decreases the total energy of the system if 
Eds + E s < 0, and we obtain the following condition of its existence 

1 /3 

T c > (h 2 adf/d s ^j = h (v/h) 1/3 . (87) 

Due to the small factor (cr/h) << 1 this condition is less restrictive than the paramagnetic limit (T c > O.66/1). 
Nevertheless the conditions of the formation of the domain structure remain rather stringent. To minimize the df/d s 
ratio (and so the effective exchange field) it is better to choose the largest possible d s thickness. However, the 
maximum thickness of the region, where superconductivity will be affected by the presence of F layer is of the order 
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of £ s . Then, even in the case of the bulk superconductor d" lax ~ £, s and the condition of the domain phase formation 
in such a case reads 

T c >h^(a/h) 1/3 . (88) 

We may conclude that for the domain phase observation it is better to choose a superconductor with a large coherence 
length £ s and the ferromagnet with low Curie temperature and small energy of the domain walls. 

The transition into the domain state is a first order one, and as all transitions related with the domain walls, it 
would be highly hysteretic. This circumstance may strongly complicate its experimental observation. To overcome 
this difficulty, it may be helpful to fabricate the S/F bilayer with a fcromagnet with a low Curie temperature O < T c . 
In such case, from the very beginning we may expect the appearence of the non-uniform magnetic structure below 
6. This system in many senses would be analogous to the ferromagnetic superconductors ErRli4B4, HoMo 6 Ss and 
HoMo 6 Se 8 . 

Bcrgcrct ct al. (2000) argued that the appearance of a nonhomogeneous magnetic order in a F film deposited on 
the bulk superconductor occurs via the second order transition and the period of the structure goes to infinity at the 
critical point. They considered the hclicoidal magnetic structure with a wave vector Q and the magnetic moment 
lying in the plane of the film. The increase of the magnetic energy due to the rotation of the moments was taken to be 
proportional to Q 2 . However, the considered magnetic structure is known to generate the magnetic field at distance 
~ Q^ 1 from the film. The contribution coming from this field makes the magnetic energy to be proportional to Q 
and not to Q 2 at a small wave- vector regime. This circumstance qualitatively change the conclusions of Bergeret et 
al. (2000) and makes the transition into a nonhomogeneous magnetic state a first-order one. 

The experiments of Miihgc et al, (1998) on the ferromagnetic resonance measurements in the Fe/Nb bilayers 
revealed some decrease of the effective magnetisation below T c for the bilayers with df < 1 nm. This thickness is 
compatible with the estimate Eq. (88), but the analysis of these experimental data by Garifullin (2002) reveals the 
possibility of the formation of islands at a small thickness of Fe layer, which may strongly complicate the interpretation 
of experimental results. 



C. Negative domain wall energy 

In the previous analysis, the energy of the domain walls was considered to be constant independent of the presence of 
the superconducting layer. It is a good approximation for a thin domain wall w << £ s - However, the phenomenon of 
superconductivity localized near the domain walls is the manifestation of the local enhancement of the superconducting 
condensation energy, which may give a negative contribution to the domain wall energy. We estimate this effect for 
a thick w >> £ s domain wall. The effect is maximum for the S/F bilayer with the relative variation of the critical 
temperature (T c —T*)T C ~ 1 at d s ^ £ s . We will suppose these conditions to be satisfied. Following the same 
reasoning as in the case of the domain wall superconductivity, we may estimate the relative local decrease of the 
pair-breaking parameter as 5 (t^ 1 ) /t^ 1 ~ (£ s / w ) 2 - Therefore the local negative contribution to the domain wall 
energy (per its unit length) coming from the superconductivity reads 

SE S ~ ~N(0)A 2 (Zs/w) 2 wd s . (89) 

The proper magnetic energy of the domain wall is Eow ~ df (a/7ra 2 ) , and for a large domain wall a ~ 8 (a/w). 
The condition of the vanishing of the total energy of the domain wall SE S + Euw — gives 

T 2 f 3 a 

-JLStL^d f a~e-df, (90) 
E F wa 3 w 3 ' K ' 

where the estimate d s ^ £ s is used. Finally, we may conclude that the energy of the domain wall may be negative for 
the system with 

Tc - e yf' (91) 

where / is the electron mean free path. We have taken into account that £ s ~ and a/£o ~ T c /Ef- If the condition 
Eq. (91) is fulfilled, the following scenario emerges. The decrease of the temperature below T* will decrease the energy 
of the domain walls, which are practically always present in a ferromagnet. The concentration of the domain walls 
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will increase and finally, when the domain wall energy will change its sign, the relatively dense domain structure will 
appear. The average distance between the domains walls in such a structure would be of the order of the domain wall 
thickness itself. Note that in the case of the small thickness of the domain wall the superconducting contribution to 
its energy is negligeable and instead of Eq. (91) we obtain the non-realistic condition T c > O (£s/0 ■ We have 

taken into account only the exchange mechanism of the interaction between magnetism and superconductivity. The 
orbital effect gives an opposite contribution to the domain wall energy, related with the out of plane magnetic field 
near the domain wall, which generates the screening currents in the superconducting layer. 

At the present time, there are no clear experimental evidences for the domain structure formation in S/F bilayers. 
The experiments of Miihge et al, (1998) on the ferromagnetic resonance measurements in the Fe/Nb bilayers revealed 
some decrease of the effective magnetization below T* for the bilayers with df < 1 nm. This thickness is compatible 
with the estimate Eq. (88), but the magnetic moment decreases continuously below T* . In addition the analysis of 
these experimental data by Garifullin (2002) reveals the possibility of the formation of islands at small thickness of 
iron layer thus reducing its magnetic stiffness. The condition Eq. (91) is apparently fulfilled in the experiments of 
Miihge et al, (1998). Therefore the decrease of the domain wall energy may be at the origin of the observed effect. 

D. Ferromagnetic film on a superconducting substrate 

Bulaevskii and Chudnovsky (2000) and Bulaevskii et al. (2002) demonstrated that the pure orbital effect could 
decrease the equilibrium domain width in the ferromagnetic film on the superconducting substrate. The ferromagnet 
with a perpendicular magnetic anisotropy is either an insulator, or it is separated from the superconductor by a thin 
insulating (e. g. oxide) layer, see Fig. 24. 

In such case the ferromagnetic film and the superconductor are coupled only by the magnetic field. It is well-known 
(Landau and Lifshitz, 1982) that the positive energy of the magnetic field favors small domains, so that the stray 
field does not spread at large distance. On the other hand, the positive domains wall energy favors a large domain 
size. The balance of these two contributions gives the equilibrium domain width Iff ~ y/wdj. In the presence of a 
superconductor, the screening currents modify the distribution of the magnetic field near the S/F interface and give 
an additional positive contribution to the energy of the magnetic field. This results in the shrinkage of the domain 
width. The energy Ed of the domain structure on the superconducting substrate reads (Bulaevskii and Chudnovsky, 
2000 and Bulaevskii et al, 2002) 



Here 7 = 1/(4ttX) and In = ?jv/(47rA) are the reduced widths of domains on a superconducting and normal substrate 
respectively, and A is the London penetration depth. The minimization of Ed over I gives the equilibrium width 
of domains. In the limit A — > oo the influence of superconductivity vanishes and I — Iff. The limit A — > 0, when 
the magnetic field does not penetrate inside the superconductor was considered by Sonin (2002). In this limit the 
shrinkage of the width of the domains is maximum and I = y/2/3lff. Then we may conclude that the influence of 
superconductivity on the domain structure is not very large and it is even less pronounced in S/F bilayer when the 
thickness of the S layer becomes smaller than the London penetration depth (Daumens and Ezzahri, 2003). 

Hclseth et al. (2002) studied the change of the Bloch domain wall structure in a ferromagnetic thin film on the 
superconducting substrate with the in-plane magnetization of the domains. It occurs that the wall experiences a small 
shrinkage, which corresponds to the increase of the energy of the domain wall. 

Recently, Dubonos et al. (2002) demonstrated experimentally the influence of the superconducting transition on 
the distribution of the magnetic domains in mesoscopic ferromagnet-superconductor structures. This finding makes 
quite plausible the observation of the effect predicted by Bulaevskii and Chudnovsky (2000) and Bulaevskii et al. 
(2002). Rearrangement of the domains normally results in the resistance change in metallic ferromagnets. In this 
context Dubonos et al. (2002) noted that domain walls' displacement due to the superconducting transition could be 
the actual mechanism of the long-range resistive proximity effects previously observed in mesoscopic Ni/Al structures 
(Petrashov et al, 1999) and Co/Al structures (Giroud et al., 1998). Note also that Aumentado and Chandrasekhar 
(2001) studied the electron transport in submicron ferromagnet (Ni) in contact with a mesoscopic superconductor (Al) 
and demonstrated that the interface resistance is very sensitive to the magnetic state of the ferromagnetic particle. 



Ed ~ 3Z H — — 



(92) 




33 



X. CONCLUSIONS 



The most striking peculiarity of the proximity effect between superconductor and ferromagnet is the damped 
oscillatory behavior of the Cooper pair wave function in ferromagnet. It results in the non- monotonous dependence 
of the critical temperature of S/F bilayers and multilayers on the F layer thickness, as well as in the formation of 
"tt"— junctions in S/F/S systems. The minimum energy of the "71""— junction is realized for the phase difference ±7r, 
and a spontaneous supercurrent may appear in a circuit containing the "7r"— junction. Two possible directions of 
the supercurrent reflect the double-degenerate ground state. In contrast to the usual junction such a state is achieved 
without external applied field. The qubit (or quantum bit) is the analog of a bit for quantum computation, describing 
by state in a two level quantum system (Nielsen and Chuang, 2000). The S/F systems open a way to create an 
environmentally decoupled (so called "quiet") qubit (Ioffe et at, 1999) on the basis of the S/F/S junction. 

The "tt" — junctions allow for a realization of the concept of the complimentary logic. In the metal-oxide supercon- 
ductor logic family the combination of the semiconducting n-p-n junctions with the complimentary p-n-p ones permits 
to significantly simplify the circuitery. The similar is possible for the Josephson junctions devices and circuits when 
the "tt" — junctions are used (Terzioglu and Beasley, 1998). The logic cells with the "tt"— junctions play a role of the 
complimentary devices to the usual Josephson logic cells. 

Recently, Ustinov and Kapluncnko (2003) proposed to use the "7r"— junction as a phase shifter in the rapid single- 
flux quantum circuits. The relatively large geometrical inductance, which is required by the single-flux quantum 
storage, may be replaced by the much smaller "7r"— junction. The advantage of the implementation of the "tt" — 
junctions is the possibility to scale the dimension of superconducting logic circuits down to the submicron size. In 
addition, the use of the "tt" — junction as a phase shifter substantially increases the parameter margins of the circuits. 

As it has been discussed in Section III.D the exchange interaction strongly affects the Andreev reflection at the 
F/S interface presenting a powerful tool to probe ferromagnets and measure their spin polarization. 

The structures consisting of "0" and "tt" — Josephson junctions can exhibit quite unusual properties. Bulaevskii et at 
(1978) demonstrated that the spontaneous Josephson vortex carrying the flux $0/2 appears at the boundary between 
"0" and "tt" —junctions. A periodic structure consisting of small (comparing with Josephson length) alternating "0" 
and "tt"— Josephson junctions may have any value of an equilibrium averaged phase difference (po in the interval 
— tt < ipo < tt, depending on the ratio of lengths of "0" and "tt"— junctions (Mints, 1998; Buzdin and Koshelev, 
2003). The S/F heterostructures provide the possibility of the realization of such "<p" — junction with very special 
two maxima current-phase relation and Josephson vortices carrying partial fluxes $0 ('A)/ 71 ") and $0 (1 — Po/n)- 

The possibility to combine in a controlled manner paramagnetic and orbital mechanisms of the interaction between 
superconductivity and magnetism makes the physics of S/F heterostructures quite rich and promising for potential 
applications. Let us mention in this context the recent observation of strong vortex pinning in S/F hybrid structures, 
the spin valve effect in F/S/F systems and the domain wall superconductivity, which open a large perspective to 
the creation of new electronics devices. The progress of controllable fabrication of high-quality heterostructures and 
especially the high-quality interfaces was crucial for the recent breakthrough in this domain. Further development of 
the microfabrication technology permits to expect another interesting findings in the near future. 
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XII. APPENDIX: 



A. Bogoliubov-de Gennes equations 

As the characteristic length of the induced superconductivity variation in a ferromagnet is small compared with a 
superconducting length, it implies the use of the microscopic theory of superconductivity to describe the proximity 
effect in S/F structures. The very convenient microscopical approach to study the superconducting properties in the 
ballistic regime (the clean limit) in the presence of spatially varying field is the use of the Bogoliubov-de Gennes 
equations (de Gennes, 1966a). The equations for electron and hole wave functions uj (r) and v± (r) are 

(H - h (r)) u T (r) + A (r) Vl (r) = (r) (93) 
A* (r) U] (r) -(H + h (r)) v L (r) = E ]V[ (r) , 

where is the quasiparticle excitation energy, Hq = —ti 2 j^ — Ep is the single particle Hamiltonian, h (r) is the 
exchange field in the ferromagnet, and the spin quantization axis is chosen along its direction. The equations for the 
wave functions with opposite spin orientation u± (r) and (r) and the excitation energy E± are obtained from Eq. 
(93) by the substitution h — -> —h. Note that the solution (ui,v-\) with energy E± may be immediately obtained from 
the solution of Eq. (93), if we choose U[ = = — uj_ and E± = —E^. The pair potential in the superconductor is 

determined by the self-consistent equation 

A(r) = A ]T « T (r) v\ (r) (1 - 2/ (£7 T )) , (94) 

Ei>0 

where / (E) is the Fermi distribution function / (E) = 1/ (1 + exp (E/T)) , and A is the BCS coupling constant. 

Assuming that the Cooper pairing is absent in the ferromagnet, we have A(r) = there. The situations when it is 
possible to obtain the analytical solutions of the Bogoliubov-de Gennes equations with spatially varying pair potential 
are very rare. However, these equations provide a good basis for the numerical calculations to treat different aspects 
of S/N and S/F proximity effects. 

B. Eilenberger and Usadel equations for ferromagnets 

Another microscopical approach in the theory of superconductivity uses the electronic Green's functions. The 
Green's functions technique for superconductors has been proposed by Gor'kov who introduced in addition to the 
normal Green's function G(r 1 ,r 2 ) the anomalous (Gor'kov) function F(r 1 ,r 2 ) (see, for example, Abrikosov et at, 
1975). This technique is a very powerful tool, but the corresponding Green's functions in a general case occur to 
be rather complicated and oscillate as a function of the relative coordinate r 1 — r 2 on the scale of the interatomic 
distance. On the other hand, the characteristic length scales for superconductivity in S/F systems are of the order 
of the layers thicknesses or damping dacay length for the induced superconductivity and, then, they are much larger 
than the atomic length. This smooth variation is described by the center of mass coordinate r = (r 1 + r 2 ) /2 in the 
Green's functions. The very convenient quasiclassical equations for the Green's functions averaged over the rapid 
oscillations on the relative coordinate has been proposed by Eilenberger (1968) (and also by Larkin and Ovchinnikov 
(1968)). 

Eilenberger equations are transport-like equations for the energy-integrated Green's functions f(r,u>, n) and 
g(r,u),n), depending on the center of mass coordinate r, Matsubara frequencies u = 7rT(2n+l) and the direc- 
tion of the unit vector n normal to the Fermi surface. For the case of S/F multilayers we may restrict ourselves to 
the situations when all quantities only depend on one coordinate x, chosen perpendicular to the layers. Introducing 
the angle 6 between the x axis and the direction of the vector n (direction of the Fermi velocity), we may write the 
Eilenberger equations in the presence of the exchange field h(x) in the form (see, for example Bulaevskii et al. (1985) 
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and a recent review on the physics of Josephson junctions by Golubov et al. (2004)) 



oj + ih(x) + — G(x, oj) ) / (x, 9, oj) + -vp cos 9^ 



^{x) + ^F{x,oj) S jg{x,9,oj), (95) 

G(x, oj) = J^-g (x, 9, oj) , F(x, oj) = J ^-f (x, 9, oj) , 

f(x,9,oj)f+ (x,9,oj)+g 2 (x,9,uj) = 1, 

where the function f + (x, n,oj) satisfies the same equation as f(x, — n,oj) with A — > A* and the presence of impurities 
is descibed by the elastic scattering time r = l/vf. The functions G(x,u>) and F{x,oj) are the Green's functions 
averaged over the Fermi surface. The Eilcnberger equations arc completed by the self-consistency equation for the 
pair potential A(x) in a superconducting layer : 

A(x) = itT\^F{x,oj). (96) 

UJ 

The BCS coupling constant A is spatially independent in a superconducting layer, while in a ferromagnetic layer it is 
equal to zero. In a superconducting layer, the self-consistency equation may also be written in the following convenient 
form 

A(x)ln^+nTj2{^--F(^))=0, (97) 

where T c0 is the bare transition temperature of the superconducting layer in the absence of proximity effect. 

Note that the presented form of the Eilenberger equations implies the natural choice of the spin quantization axis 
along the direction of the exchange field, and the only difference with the standard form of these equations is the 
substitution of the Matsubara frequency wbyw + ih(x). 

Usually, the electron scattering mean free path in S/F/S systems is rather small. In such a dirty limit, the angular 
dependence of the Green's functions is weak, and the Eilenberger equations can be replaced by the much simpler 
Usadel (1970) equations. In fact, the conditions of the applicability of the Usadel equations are T c t <C 1 and hr <C 1. 
The second condition is much more restrictive due to a large value of the exchange field (h 3> T c ). The Usadel 
equations only deal with the Green's functions G(x, oj) and F(x,oj) averaged over the Fermi surface : 



d 2 d 2 
G(x,oj, h)-^F (x,u, h) - F(x,oj, h)-^G {x,uj, h) 



D 

+ (lu + ih{x)) F(x, oj, h) = A(x)G (x, oj, h) , (98) 



G 2 (x,oj,h) + F(x,oj,h)F + (x,h,oj) = 1, 

D = ^vpl is the diffusion coefficient which is different in S and F regions and the equation for the function F + (x, h, oj) 
is the same as for F{x,oj, h) with the substitution A — > A*. Here also the only difference with the standard form of 
the Usadel equations is the substitution oj by oj + ih(x). 

The equations for the Green's functions in F and S regions must be completed by the corresponding boundaries 
conditions at the interfaces. For the Eilenberger equations they were derived by Zaitsev (1984) and for the Usadel 
equations by Kupriyanov and Lukichev (1988). These boundary conditions take into account the finite transparency 
(resistance) of the interfaces - see Eq. (32). 

The most important pair-breaking mechanism in the ferromagnet is the exchange field h. However a disorder in 
the lattice of magnetic atoms creates centers of magnetic scattering. In ferromagnetic alloys, used as the F layer in 
S/F/S Josephson junctions, the role of the magnetic scattering may be quite important. Note that even in the case of 
a perfect ordering of the magnetic atoms, the spin-waves will generate magnetic scattering. The natural choice of the 
spin-quantization axis used implicitly above is along the direction of the exchange field. The magnetic scattering and 
spin-orbit scattering mix up the up and down spin states. Therefore to describe this situation it is needed to introduce 
two normal Green's functions G\ <~ (ip^ipt), G2 ~ (ip±ipt) and two anomalous ones F\ ~ (tp^tpi), F 2 <~ (V'tV't) • 
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The microscopical Green's function theory of superconductors with magnetic impurities and spin-orbit scattering was 
proposed by Abrikosov and Gorkov (1960, 1962). The generalization of the Usadel equations (98) to this case gives 



D 
~2 



d 2 d 2 
G\-^—~F\ — Fi—^rGi 
ox ox 



+ [ UJ + th+[— + —)G 1 )F 1 + 



G 1 (F 2 - Fx) (— - —) + F 1 (G 2 - Gi) (- + —)= A(x)G u 

G\ {x,lj, h) + Fi(x,w, h)F^{x, h,uj) = 1, 

and the similar equation for F 2 with the indices substitution 1 <-» 2. Here is the spin-orbit scattering rate, while 
the magnetic scattering rates are t~ x = t 2 x (S 2 ) / S 2 and t^T 1 = t 2 x (S 2 ) /S 2 . The rate t 2 x describes the intensity 
of the magnetic scattering via exchange interaction and we follow the notation of the paper of Fulde and Maki (1966). 
In the spatially uniform case the equations (??) are equivalent to those of the Abrikosov- Gorkov theory (1960, 1962) 
(see also Fulde and Maki, 1966). Demler et al. (1997) analyzed the influence of the spin-orbit scattering on the critical 
temperature of the S/F multilayers. The equations (Demler et al, 1997) corresponds to the limit A — > 0, Gi 2 = 1 in 
(??). 

The ferromagnets used as F layers in S/F heterostructures reveal strong uniaxial anisotropy. Then the magnetic 
scattering in the plane (xy) perpendicular to the anisotropy axis is negligeable. Moreover due to the relatively small 
atomic numbers of the F layers atoms the spin-orbit scattering is expected to be weak. In such case there is no 
spin mixing scattering anymore and the Usadel equations retrieve the initial form (98) with the substitution of the 
Matsubara frequencies by lu — > lu + G/t s , where r s _1 = t~ x = t^ 1 (S 2 ) / S 2 may be considered as a phenomenological 
parameter describing the intensity of the magnetic scattering (Buzdin, 1985). 

The linearized Usadel equation in the ferromagnet reads 

1\ Dfd 2 F f 

u \+ih8gn(u) + -)F f --+-^£=0. (99) 

If t s T c « 1, we may neglect \u\ in Eq. (99) and the exponentially decaying solution has the form 

Ff (x,u > 0) = Aexp(-x{h +ik 2 )), (100) 



with k\ = -^j v \7T+ a 2 + a and k 2 = y \7T+ a 2 — a, where a = l/(r s /i). In the absence of magnetic scattering, 
the decaying and oscillating wave vectors are the same k\ = k 2 . The magnetic scattering decreases the characteristic 
decaying length and increases the period of oscillations. In practice, it means that the decrease of the critical current 
of S/F/S junction with the increase of df will be more strong. Note that the spin-orbit scattering (in contrast to the 
magnetic scattering) decreases the pair-breaking effect of the exchange field (Demler et al., 1997) and both scattering 
mechanisms decrease the amplitude of the oscillations of the Cooper pair wave function. In some sense the spin-orbit 
scattering is more harmful for these oscillations because they completely disappear at r^ 1 > h. The observation on 
experiment of the oscillatory behavior of T c in S/F multilayers is an indirect proof of the weakness of the spin-orbit 
scattering. 

The expression for I c (2df) dependence (54) may be generalized to take into account the magnetic scattering 



I c R n = 64 — Re 



,7rT 

e 



E 



2g a ,yexp(-2g a ,2/)$ tl 



[va-^c + i + i 



(101) 



where the functions 



(O + w)- 



-, = ^/2i + 2a + 2u/h, rfc 



a 



a + i + u/h 



(102) 



Near T c and in the limit h » T c and 2dfk 2 >> 1 it possible to obtain the following simple analytical expression 
for the critical current 
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7r5cr/A 2 fci 
2eT c 



k 

cos (2djk 2 ) + — sin (2dfk 2 ) 
k\ 



exp(-2d/fci). (103) 



We see that due to the magnetic scattering the decaying length of the critical current £fi — l/k\ may be substantially 
smaller than the oscillating length £/2 = l/fa- 

As it has been noted above, the condition of the applicability of the Usadel equations, hr <C 1, is rather restrictive 
in ferromagnets due to the large value of the exchange field. Therefore, it is of interest to retain in the Usadel 
equations the first correction in the parameter hr. The first attempts to calculate this correction were made by 
Tagirov (1998) and Proshin and Khusainov (1998) and resulted in the renormalization of the diffusion constant of 
the F layer Df — > Df(l — lihrsign (to)). Later on, the similar renormalization has been proposed by Bergeret et 
al. (2001c) and Baladie and Buzdin (2001). The critical analysis of this renormalization by Fominov et al. (2002) 
(see also Fominov et al. 2003b and Khusainov and Proshin, 2003) revealed the inaccuracy of this renormalization, 
but did not provided the answer. The careful derivation of the Usadel equation for an F layer retaining the linear 
correction over the parameter hr was made by Buzdin and Baladie (2003) and simply resulted in a somewhat different 
renormalization of the diffusion constant Df — > Df(l — OAihrsign (uj)). The coefficient in the parameter hr occurs 
to be rather small which provides more confidence in the description of F layers in the framework of the Usadel 
equations. Note that this renormalization of the diffusion constant increases the decaying characteristic length and 
decreases the period of oscillations, which is opposite to the influence of the magnetic scattering. 

The Usadel equations give the description of Green's functions only on average. Zyuzin et al. (2003) pointed out 
that due to the mcsoscopic fluctuations, the decay of the anomalous Green's function Ff at distances much larger 
than £f is not exponential. In result, the Josephson effect in S/F/S systems may be observed even with a thick 
ferromagnetic layer. 

The Eilenberger and Usadel equations adequately describe the weak ferromagnets, where h « Ep and the spin-up 
«FT, an d spin-down vf\, Fermi velocities are the same. When the parameters of the electrons spectra of the spin- 
up and spin-down bands are very different, the quasiclassical approach fails. However, if the characteristics of the 
spin bands are similar, the Eilenberger and Usadel equations are still applicable. Performing the derivation of the 
Eilenberger equation in such case, it may be demonstrated that the Fermi velocity vf in Eq. (95) must be substituted 
by [vf] + v fi) /2 and the scattering rate 1/t by (1/tj + 1/tj.) /2. In consequence, the diffusion coefficient Df in the 
Usadel equation becomes (1/6) (vp^ + vfi) 2 / (1/ T T + 1/ T l) ■ Let us stress that such renormalization is justified only 
for close values of VF^&nd Vf\, (as well as T| and rjj. Otherwise the Bogoliubov-de Gennes equations must be used 
for the description of the proximity effect in the strong ferromagnets. 
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TABLE I. Characteristic length scales of S/F proximity effect. 



Thermal diffusion length Lt 


V 2¥T 


Superconducting coherence length £ s 


2 F r£, in pure limit 
\j ~TkT~ m dirty limit 


Superconducting correlations decaying length £1/ in a ferromagnct 


H£ in pure limit 

if = \J ^ in dirty limit 


Superconducting correlations oscillating length £2/ in a ferromagnet 


^J- in pure limit 

S,f = sf^f- m dirty limit 



Figure captions 



FIG. 1. The (T,H) phase diagram for 3D superconductor. At temperature below T* = 0.56T C the second order 
transition occurs from the normal to the non-uniform superconducting FFLO phase. The dashed line corresponds 
to the first order transition into the uniform superconducting state, and the dotted line presents the second order 
transition into the uniform superconducting state. 

FIG. 2. Energy band of ID superconductor near the Fermi energy. Due to the Zeeman splitting the energy of the 
electrons with spin orientation along the magnetic field (|) decreases - dotted line, while the energy of the electrons 
with the opposite spin orientation (J.) increases - dotted line. The splitting of the Fermi momenta is ±5kF, where 
SttF = hbH/vf ■ The Cooper pair comprises one electron with the spin (|) and momentum k F + 5k Fl and another 
electron with the spin ({) and momentum — k F + Sk F . The resulting momentum of the Cooper pair is non-zero : 
k F + 5k F + {-k F + 5k F ) = 2Sk F ^ 0. 

FIG. 3. Schematic behavior of the superconducting order parameter near the interface (a) superconductor-normal 
metal, and (b) superconductor- ferromagnet. The continuity of the order parameter at the interface implies the absence 
of the potential barrier. In general case at the interface the jump of the superconducting order parameter occurs. 

FIG. 4. Measurements of the differential conductance by Kontos et al. (2001) for two Al/Al 2 03/PdNi/Nb junctions 
with two different thicknesses (50 A and 75 A) of the ferromagnetic PdNi layer. A 1500-A-thick aluminium layer was 

evaporated on SiO and then quickly oxidized to produce a AI2O3 tunnel barrier. Tunnel junction areas were defined 
by evaporating 500 A of SiO through masks. A PdNi thin layer was deposited and then backed by a Nb layer. 

FIG. 5. Experimental data of Jiang et al. (1995) on the oscillation of the critical temperature of Nb/Gd multilayers 
vs thickness of Gd layer da for two different thicknesses of Nb layers : (a) d^b —600 A and (b) d^b =500 A. Dashed 
line in (a) is a fit by the theory of Radovic et al. (1991). 

FIG. 6. S/F multilayer. The axe x is chosen perpendicular to the planes of S and F layers with the thicknesses 2d s 
and 2df respectively, (a) The curve ^(2) represents schematically the behavior of the Cooper pair wave function in 
"0"- phase. Due to the symmetry reasons the derivative of 'J (and F) is zero at the centers of S and F layers. The 
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case of the "0"- phase is equivalent to the S/F bilayer with S and F layers thicknesses d s and df respectively, (b) 
The Cooper pair wave function in "7r"- phase vanishes at the centers of F layers and ^(x) is antisymmetric toward 
the center of F layer. 

FIG. 7. The dependence of the critical temperature on the thickness of F layer for "0" -phase (solid line) and "tt"- 
phase (dotted line) in the case of the transparent S/F interface. Note that the highest transition temperature T* 
corresponds to the lowest point. The dimensionless thickness of F layer 2y — 2df /£/ and the first transition from 
"0"- to "7r"- phase occurs at 2c?/ = 2. 36£/. The parameter To = —jj^fj- 

FIG. 8. The critical temperature of "0"- phase (solid line) and 'V- phase (dashed line) as a function of the 
dimensionless thickness of F layer 2y = 2d//£/ for different S/F interface barriers 7 = 75 • 

(a) The dimensionless pair-breaking parameter fo = 47rT c 2d ^ f ga. =21. 

(b) The dimensionless pair-breaking parameter fo = 20.05. 

FIG. 9. Variation of the critical temperature of the Nb/Cuo. 43N10.57 bilayer with the F layer thickness (Ryazanov 
et al. 2003). Theoretical fit (Fominov et al., 2002) gives the exchange field value h ~ 130 K and the interface 
transparency parameter 73 ~ 0.3. 

FIG. 10. Geometry of the S/F/S junction. The thickness of the ferromagnetic layers is 2df and the both S/F 
interfaces have the same transparencies, characterized by the coefficient 73. 

FIG. 11. Critical current of the S/F/S Josephson junction near T c as a function of the dimensionless thickness of 
F layer 2y = 2d//£/. There are no barriers at the S/F interfaces (73 = 0), R n is the resistance of the junction and 

T 7 _ 7rA 2 

FIG. 12. Temperature dependences of the critical thickness 2dj of F layer, corresponding to the crossover from 
"0"- to "71""- phase in the limit of very small boundary transparency for different values of the exchange field. 

FIG. 13. Non-monotonous temperature dependences of the normalized critical current for low transparency limit: 
curve 1: h/T c = 10 and 2d//£/ = 0.84; curve 2: h/T c = 40 and 2d//£/ = 0.5; curve 3: h/T c = 100 and 2d//£/ = 0.43. 

FIG. 14. Critical current I c as a function of temperature for Cuo.4sNio.52 junctions with different F layers thicknesses 
2dp. At the thickness of F layer of 27 nm the temperature mediated transition between "0"- and 'V- phases occurs. 
Adapted from (Ryazanov et al., 2001a). 

FIG. 15. Critical current I c at T = 4.2 K of Cuo.47Nio.53 junctions as a function of F layer thickness (Ryazanov et 
al, 2005). Two "0" — 'V transitions are revealed. The theoretical fit corresponds to the Eq. (101) in Appendix B, 
taking into account the presence of the magnetic scattering with parameters a = l/(r s h) = 1.33 and £j = 2.4 nm. 
The inset shows the temperature mediated "0" — "71"" transition for the F layer thickness 11 nm. 

FIG. 16. The experimental points correspond to the measurements of the critical current, done by Kontos et 
al. (2002) vs the PdNi layer thickness. The theoretical curve is the fit of Buzdin and Baladie (2003). The fitting 
parameters are ~30 A and ~ 110 /j,V. 

FIG. 17. Experiments of Guichard et al. (2003) on the diffraction pattern of SQUID with "0"- and "7r" -junctions. 
There is no shift of the pattern between a " — 0" and " tt — 7r" SQUIDs. The <&o/2 shift is observed between a " — 7r" 
and "0 — 0" or "tt — 7r" SQUIDs. The "0"- and 'V -junctions were obtained by varying the PdNi layer thickness. 

FIG. 18. Earlier observation by Deutscher and Meunier (1969) of the spin- waive effect on In film between oxidized 
FeNi and Ni layers. The figure presents the resistive measurements of the critical temperature in zero field: dashed 
line, after application of 1 T field parallel to the ferromagnetic layers; solid line, after application of the -1 T field and 
subsequently +0.03 T field to return the magnetization of FeNi layer. 

FIG. 19. Geometry of the F/S/F sandwich. The thickness of S layer is 2d s and two F layers have identical 
thicknesses df. 

FIG. 20. Influence of the S/F interface transparency (parameter 7 = jb (£«/£/)) on tnc ^c* vs df dependence 
(Baladie and Buzdin, 2003). The thickness of F layer is normalized to the £/. The dimensionless pair-breaking 
parameter f = AnT c 2d jf } is chosen constant and equal to 4. The full line corresponds to the antiparallel case, and 
the dashed line to the parallel case. One can distinguish four characteristic types behavior: (a) weakly non-monotonous 
decay to a finite value of T* , (b) reentrant behavior for the parallel orientation, and (c) and (d) monotonous decay 
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to T* = with (d) or without (c) switching to a first-order transition in the parallel case. In (d), the dotted line 
presents schematically the first order transition line. 

FIG. 21. The calculate dependence of the superconducting transition temperature vs inverse reduced half-thickness 
d*/d s of the superconducting layer for parallel and antiparallcl alignments for the transparent interface (75 = 0) and 

1 /2 

thick ferromagnetic layer (df » £/). The effective length d* = (<Tf/a s ) (D s /4ttT c ) (h/Df) 1 . 

FIG. 22. The (T, h) - phase diagram of the atomic S/F multilayer in the limit of the small transfer integral t «T C . 

FIG. 23. S/F bilayer with domain structure in the ferromagnetic layer. The period D of the domain structure 
(D = 2tt/Q) is smaller than the superconducting coherence length £ s . 

FIG. 24. The ferromagnetic film with perpendicular anisotropy on a superconducting substrate. 
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